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ABSTRACT 


During  the  3  years  of  this  ARO  Grant  we;  Measured  and  coMputed  the  drop 
size  in  the  imnediate  vicinity  of  the  nozzle  in  Diesel-type  sprays  (i'8); 
nea^red  (Appendix  B  and  Ref.  9)  and  computed  (Appendix  C  and  Ref.  10)  the  drop 
velocities  in  the  farfield  of  non  vaporizing  Diesel-type  sprays;  Measured 
(and  are\;onputing)  the  drop  velocities  in  the  farfield  of  vaporizing  Diesel- 
type  sprays.  The  everal-^ conclusions  are  that:  the  breakup  of  full-cone  Diesel 
Jets  is  due  primarily  to  aerodynamic  interaction  between  the  liquid  surface  and 
the  chamber  gas;  the  initial  average  size  and  velocity  of  the  drops  formed  by 
the  breakup  of  the  outer  surface  of  the  Jet  can  now  be  estimated  and  so  can  the 
length  of  the  Intact  core  but  with  greater  uncertainty;  after  breakup  a  steady 
full-cone  spray  undergoes  a  development  similar  to  that  of  Incompressible  Jets; 
in  Diesel  engines  it  is  the  development  region  that  is  of  Importance;  the  major 
parameter  in  the  difference  between  full-cone  sprays  and  incompressible  Jets  is 
the  liquid- to-gas  density  ratio;  in  the  absence  of  vaporization  the  most 
advanced  of  the  available  spray  model  reproduce  full-cone  sprays  with  adequate 
realism  for  applications;  it  is  not  yet  known  how  vaporization  changes  the 
structure  of  these  sprays  and  if  available  models  adequately  reproduce  the 


changes,  j^e  also  produced  a  number  of  summary  papers  on  the  mechanism  of  atomi¬ 
zation  (Aj^endlx  A  and  Refs.  12,13)  and  on  the  structure  of  Diesel  sprays 


(Appendix  D  and  Ref.  14).  One  student  received  his  Ph.D.  on  this  subject  and  is 


now  at  the  Research  Laboratories  of  GM  (Dr.  K.-J.  Wu)  and  several  staff  members 
and  visitors  work^  with  us  and  have  now  moved  on  to  other  positions  (Drs.  R.D. 


Reitz,  D.A.  Santavicca,  A.  Coghe,  Y.  Onuma,  and  P.  Felton)] 
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REPORT 

Pour  papers  are  given  in  the  appendices.  Each  of  the  four  papers  is  a  sum- 
■ary  of  much  work  conducted  under  this  Grant  and  the  following  paragraphs  are  a 
summary  of  the  summaries. 

In  1972-1976,  we  made  some  multi-dimensional  computations  of  engine  sprays 
(e.g.,  CST,  Vol.  8.  pp.  69-84,  1973;  SAE  TRANS.,  Vol .  84,  pp.  3317-3340;  1975; 
CST,  Vol.  12,  pp.  63-74,  1976;  AIAA  J. ,  Vol.  16,  pp.  1053-1061,  1978).  Those 
were  the  first  multi-dimensional  two-phase  engine  flow  computations  and  from 
them  we  learned  that  quantitative  knowledge  of  the  mechanism  of  break-up  of  fuel 
Jets  and  of  the  structure  of  dense  sprays  was  so  poor,  and  their  effects  so 
important,  that  the  reliability  of  even  the  computed  trends  was  questionable. 

Vast  qualitative  knowledge  of  sprays  was  accumulated  in  the  30's,  primarily 
at  NACA,  often  with  astute  experimental  techniques.  But  when  space-resolved 
computations  of  in-cylinder  events  became  necessary  and  feasible,  the 
corresponding  space-resolved  quantitative  spray  information  simply  was  missing. 

Thus,  in  1976  we  started  a  program  of  detailed  measurements  of  the  break-up 
of  fuel  jets  to  establish  its  mechanism  with  the  final  objective  of  determining 
equations  for  the  initial  drop  sizes  and  velocities.  ARO  provided  its  first 
3-year  support  for  the  period  July  78  -  June  81.  Full-cone  sprays  from  single 
cylindrical  orifices  were  injected  into  compressed  gases  in  a  constant-volume 
vessel  and  the  spray  angle  and  the  tip  penetration  velocity  were  measured  with 
photographic  techniques  in  a  large  variety  of  conditions.  Room  temperature  was 
employed  and  vaporization  was  negligible. 


During  the  second  3-year  ARO  grant  for  the  period  July  81-June  84  drop 
sizes  were  measured  by  photography  in  the  immediate  vicinity  of  the  nozzle 


exit.  Also  Measured  by  LOV  were  drop  velocities  within  nonvaporizing  sprays  at 
distances  between  2.3  cn  and  10  cm  froM  the  nozzle  (Appendix  B) .  The  diaweter 
of  the  nozzles  used  was  between  76  and  335  mm.  Although  drop  sizes  and  veloci¬ 
ties  had  been  measured  before,  they  had  not  been  measured  respectively  in  the 
Immediate  vicinity  of  the  nozzle  exit  and  within  very  dense  sprays  as  we  did. 

With  the  help  of  the  new  data,  we  sorted  out  the  large  amount  of  often- 
contradicting  hypotheses  that  had  been  made  about  the  formation  of  Diesel-type 
sprays  and  concluded  that  atomizing  liquid  jets  are  broken  up  by  aerodynamic 
forces  that  are  due  to  the  llquld-gas  relative  velocity.  These  forces  selec¬ 
tively  amplify  liquid  perturbations  that  are  generated  within  the  nozzle.  This 
particular  mechanism  was  first  proposed,  in  a  quantitative  way,  by  W.E.  Ranz 
<2).  Vie  now  understand,  or  we  think  we  do,  the  relationship  between  the  many 
modes  of  breakup  of  liquid  Jets.  They  are  different  manifestations  of  the  same 
forces  that  align  themselves  differently  under  different  conditions.  These  for¬ 
ces  act  on  perturbations  of  the  llquld-gas  Interface.  The  perturbations.  In 
turn,  are  often  generated  within  the  nozzle.  At  low  jet  velocities,  the  nozzle 
perturbations  actual-ly  are  not  very  Important,  but  they  become  progressively 
more  significant  as  the  jet  velocity  increases  (and  the  intact  length  decreased) 
and  are  essential  components  of  the  breakup  in  the  atomization  regime,  l.e.  for 
engine  fuel  jets.  The  mechanisms  of  breakup  of  round  liquid  Jets  are  discussed 
In  Appendix  A. 

Actually,  as  far  as  the  formation  of  engine  sprays  is  concerned,  we  get 
several  lucky  breaks.  Practical  nozzle  geometries  do  not  vary  much,  due  to 
construction  considerations,  and  Influence  mostly  the  initial  angle  of  the  spray 
and,  possibly,  the  solid  core  length  Both  angle  and  solid  core  length  are 


Influenced  also,  and  auch  nore  Markedly  and  Importantly ,  by  the  gas-liquid- 
density  ratio.  Significantly,  they  are  not  strongly  affected  by  any  of  the 
other  many  parameters  (Eqs.  11  and  10a  of  App.  A,  also  Eqs.  15  and  17  of  App. 

D)  except  for  the  solid  core  length  that  is  proportinal  also  to  the  nozzle 
diameter  (Eq.  17  of  App.  D) .  The  Initial  size  of  the  drops  that  are  formed  by 
the  breakup  of  the  outer  surface  of  the  liquid  jet  is  determined  mostly  by  the 
Weber  number  (Eq.  12  of  App.  A,  also  Eq.  16  of  App.  D).  The  predictions  of  the 
initial  angle  has  now  been  supported  by  numerous  measurements.  That  of  the  ini¬ 
tial  drop  size  has  been  evaluated  only  preliminarily  But  the  prediction  of 

the  Intact  core  length  is  both  weaker  and  largely  untested.  It  is  weaker 
because  it  is  based  on  the  questionable  assumption  that  the  stability  con¬ 
siderations  that  apply  to  the  outer  surface  of  the  liquid  also  apply  unchanged 
within  the  Jet  and  all  the  way  down  to  its  axis.  It  is  largely  untested  because 
the  length  of  the  liquid  core  has  been  measured  only  once  and  with  a  tech¬ 
nique  that  can  be  criticized. 

Getting  back  to  our  overview  and  assuming  that  we  can  estimate  the  initial 
angle  of  the  spray,  the  initial  size  of  the  drops,  and  the  shape  and  length  of 
the  solid  core,  l.e.  the  drop  generating  surface,  we  can  then  compute  the  sub¬ 
sequent  development  of  the  spray.  In  particular,  we  can  compute  the  farfield 
quantities  that  we  measured  (Appendix  B),  that  is:  the  axial  and  radial  velocity 
distributions  of  the  drops  at  all  radial  and  axial  positions  within  the  spray 
and  for  a  variety  of  liquid  and  gas  conditions. 

The  computations  and  comparisons  with  the  measured  drop  velocity  are 
discussed  in  Appendix  C.  They  clearly  show  that  our  model  for  Diesel  sprays 


accurately  predicts  the  average  values  of  the  sprays  parameters  but  not  the  the 


fluctuating  conponent  of  the  drop  velocity.  More  work  Is  now  needed  on  the 
interaction  between  droplet  notion  and  gas  turbulence.  In  Appendix  D,  the 
fourth  appendix,  we  have  sunmarized  what  we  have  learned  about  the  fornation  and 
structure  of  full-cone  sprays  in  all  of  these  studies. 

All  the  neasureaents  and  conputations  that  we  have  discussed  so  far  were 
■ade  under  Diesel-type  conditions  except  for  the  temperature  that  had  its  roon 
value.  It  has  been  advisable  to  do  so  to  avoid  the  trap  of  having  too  nany  pro¬ 
cesses  to  sort  out  sinultaneously .  Actually  we  have  now  conpleted  a  second 
round  of  drop  velocity  neasureaents  similar  to  those  discussed  in  Appendix  B 
but  in  vaporizing  sprays.  We  are  in  the  process  of  studying  this  new  experimen¬ 
tal  information  and  then  we  intend  to  make  a  second  set  of  conputations  and  com¬ 
parisons  similar  to  those  of  Appendix  C. 

Notice,  however,  that  to  nedce  measurements  in  vaporizing  sprays,  first  we 
had  to  modify  our  current  experimental  apparatus.  It  is  difficult  to  obtain 
controlled,  repeatable  gas  and  liquid  environments  above  room  temperature  with  a 
closed  gas  vessel  of  limited  size.  Vaporization  and  wall  heat  losses  cool  the 
gas  and  condensation  on  windows  limits  optical  access,  at  very  least  compli¬ 
cating  repeatability.  Also,  the  need  to  empty  the  vessel  after  each  injection 
conflicts  with  the  long  time  required  to  gather  the  large  amount  of  data  asso¬ 
ciated  with  detailed  measurements.  Thus  we  decided  to  build  a  continuous  flow 
system  and,  in  order  to  control  separately  liquid  and  gas  temperature,  it  was 
also  necessary  to  redesign  the  fuel  injection  system  to  allow  for  its  water 
cooling. 

Much  work  was  al  j  done  on,  and  with,  the  spray  model  to  which  we  have 
often  referred.  This  model  is  important  because  it  represents,  i.e.  embodies, 


our  current  understanding  of  the  structure  of  these  sprays  and  because  Is 
already  In  use  In  several  industries,  including  Cummins  and  Komatsu. 

The  model  computes  the  unsteady  penetration  and  the  steady  state  of  axlsym- 
metrlc  dense  sprays  using  information  about  drop  sizes  and  velocities  from  the 
break-up  of  the  fuel  jet  at  the  nozzle  exit  as  input.  The  model  solves  the 
coupled  conservation  equations  for  the  gas  and  liquid  with  an  Eulerlan  scheme 
for  the  gas  and  a  Lagranglan  scheme  for  the  drops.  The  model  includes  exchange 
functions  for  droplet  coalescence,  grazing  collisions  and  accounts  for  the 
effect  of  liquid  volume  fraction  on  the  rates  of  exchange  of  mass,  momentum  and 
energy  between  the  gas  drops.  The  functions  representing  these  thick  spray  phe¬ 
nomena  have  been  deduced  from  experimental  and  theoretical  work  in  cloud  phy¬ 
sics,  fluidized  and  packed  beds,  and  transport  processes  in  fixed  arrays  of 
spheres . 

In  the  model,  the  gas  flow-field  and  droplet  size,  velocity  and  temperature 
probability  distribution  functions  are  computed  as  functions  of  position  and 
time  in  the  spray.  The  optically  thick  region  of  the  spray  close  to  the  nozzle 
exit  is  also  computed.  This  region  of  dense  sprays  is  presently  Inaccessible  to 
standard  probing  measurements  and  the  model  offers  a  means  to  bridge  the 
measurements  of  the  jet  break-up  process  at  the  nozzle  exit  to  the  LDV  measure¬ 
ments  of  droplet  velocity  and  size  distributions  within  the  optical ly-thinner 
downstream  spray. 

At  present  there  are  three  advanced  spray  models;  the  CONCHAS-SPRAY  of  the 
Los  Alamos  National  Laboratory,  the  RPM-SPRAY  of  Gosmcn,  and  ours.  These  models 
differ  both  in  their  physics  and  in  their  numerics.  But  they  also  overlap. 

Thus  Gosman  and  we  use  k-e  gas  turbulence  submodels  whereas  LANL  uses  a  k-sub- 


grid.  For  full-cone  sprays,  we  all  use  the  initial  conditions  from  the  Jet 
break-up  process  determined  at  Princeton.  We  Include  computations  of  drop 
collisions  and  coalescence  and  of  the  effect  of  the  gas-volume  fraction  on  the 
gas-liquid  rates  of  exchange  of  mass,  momentum,  and  energy.  Gosman  does  not 
consider  such  processes  and  LANL  has  recently  adopted  our  collisions  and 
coalescence  submodel.  We  all  use  the  stochastic  technique  of  LANL  for  the  com¬ 
putation  of  drop  events,  but  somewhat  different  numerical  techniques  for 
the  solution  of  the  equations. 

We  believe  that  this  diversity  is  very  healthy.  It  allows  for  the  eva¬ 
luation  of  different  approaches  while  tested  best  features  are  eventually 
adopted  by  all.  Moreover  the  three  efforts  complement  each  other  very  well 
through  differences  in  emphasis.  We  have  concentrated  on  detailed  testing  of 
the  physics  in  our  model  at  the  expense  of  model  flexibility.  LANL  has  adopted 
the  greatest  flexibility  and  its  CONCHAS-SPRAY  is  the  most  general  of  the  codes 
for  engine  computations.  Gasman's  model  falls  in  between  CONCHAS-SPRAY  and  our 
model . 

In  summary  then,  during  the  3  years  of  this  ARO  Grant  we:  measured  and  com¬ 
puted  the  drop  size  in  the  immediate  vicinity  of  the  nozzle  in  Diesel-type 
sprays  (I-®);  measured  (Appendix  B  and  Ref.  9)  and  computed  (Appendix  C  and  Ref. 
10)  the  drop  velocities  in  the  farfield  of  non  vaporizing  Diesel-type  sprays; 
measured  (and  are  computing)  the  drop  velocities  in  the  farfield  of 

vaporizing  Diesel-type  sprays.  The  overall  conclusions  are  that:  the  breakup  of 
full-cone  Diesel  jets  is  due  primarily  to  aerodynamic  interaction  between  the 
liquid  surface  and  the  chamber  gas;  the  initial  average  size  and  velocity  of  the 
drops  formed  by  the  breakup  of  the  outer  surface  of  the  jet  can  now  be  estimated 
and  so  can  the  length  of  the  intact  core  but  with  greater  uncertainty;  after 
breakup  a  steady  full-cone  spray  undergoes  a  development  similar*  to  that  of 


incompressible  Jets:  in  Diesel  engines  it  is  the  development  region  that  is  of 
importance;  the  major  parameter  in  the  difference  between  full-cone  sprays  and 
Incompressible  jets  is  the  liquld-to-gas  density  ratio;  in  the  absence  of 
vaporization  the  most  advanced  of  the  available  spray  model  reproduce  full-cone 
sprays  with  adequate  realism  for  applications:  it  is  not  yet  known  how  vaporiza¬ 
tion  changes  the  structure  of  these  sprays  and  if  available  models  adequately 
reproduce  the  changes.  We  also  produced  a  number  of  summary  papers  on  the 
mechanism  of  atomization  (Appendix  A  and  Refs.  12,13)  and  on  the  structure  of 
Diesel  sprays  (Appendix  D  and  Ref.  14).  One  student  received  his  Ph.D.  on  this 
subject  and  is  now  at  the  Research  Laboratories  of  GM  (Dr.  K.-J.  Wu)  and  several 
staff  members  and  visitors  worked  with  us  and  have  now  moved  on  to  other  posi¬ 
tions  (Drs.  R.D.  Reitz,  D.A.  Santavlcca,  A.  Coghe,  Y.  Onuma,  and  P.  Felton). 
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ABSTRACT 


The  Bcchanisos  of  breakup  of  liquid  jets  injected  fr«n  a  single  hole 
orifice  into  a  gaseous  medium  are  reviewed.  As  the  jet  velocity  is 
increased,  or  as  other  operating  conditions  are  changed  appropriately, 
four  breakup  regimes  are  identified.  The  four  regimes  correspond  to 
different  combinations  of  liquid  inertia,  surface  tension  and 
aerodynamic  forces  acting  on  the  jet.  They  have  been  called  the 
Rayleigh,  the  first  wind-induced,  the  second  wind-induced  and  the 
atomization  regimes.  In  each  of  the  four  regimes,  the  outcome  is  also 
influenced  by  the  initial  state  of  the  jet.  This  influence  appears  to 
grow  in  iaiportance  with  increasing  jet  velocity.  The  existence  of  the 
regimes  is  consistent  with  a  stability  analysis  of  the  liquid  surface. 
It  appears  that  the  stability  analysis  can  account  also  for  the 
influence  of  the  initial  state  of  the  jet  but  little  work  has  been  done 
in  this  area  and  the  breakup  process  at  high  jet  velocities  is 
understood  more  qualitatively  than  quantitatively  at  present. 


1.  IHTROPUCTION 


The  mechaniSDS  of  jet  breakup  discussed  here  result  from  the  steady 
injection  of  a  liquid  through  a  single  hole  nozzle  into  a  quiescent  gas. 
The  complexity  of  the  breakup  process  is  due  to  the  unusually  large 
number  of  parameters  which  influence  It,  including  the  details  of  the 
design  of  the  nozzle,  the  jet's  velocity  and  turbulence,  and  the 
physical  and  thermodynamic  states  of  both  liquid  and  gas.  In  this  work 
we  discuss  a  framework  by  means  of  which  some  of  the  underlying 
mechanisms  of  breakup  can  be  organized  and,  eventually,  be  better 
understood. 

The  approach  followed  is  to  divide  the  jet  breakup  phenosena  of  interest 
into  various  breakup  regimes.  These  regimes  reflect  differences  in  the 
appearance  of  the  jet  as  the  operating  conditions  are  changed.  We  then 
attempt  to  relate  these  regimes  to  limiting  cases  of  a  stability 
analysis  of  liquid  jets.  The  analysis  considers  the  growth  of  initial 
perturbations  of  the  liquid  surface  and  Includes  the  effects  of  liquid 
inertia,  surface  tension,  viscous  and  aerodynamic  forces  on  the  jet.  The 
theory  is  found  to  offer  a  reasonably  coiqplete  description  of  the 
breakup  mechanisms  of  low  speed  jets.  For  high  speed  liquid  jets 
however,  the  initial  state  of  the  jet  appears  to  be  progressively  more 
important  and  less  understood.  We  will  summarize  the  results  of  recent 
research  aimed  at  closing  this  gap. 

Zt  should  be  pointed  out  that  we  consider  only  fluid-dynamic 
instabilities  and  that  there  are  other  causes  of  breakup  such  as 
superheating*,  electrostatic  charge*,  acoustical  excitation*,  and 
chemical  reactions.  Even  in  the  field  of  fluid-dynamic  instabilities 
there  are  difficulties;  many  parameters  Influence  the  outcome;  the 
magnitude  of  these  parameters  vary  over  broad  ranges;  there  are  no 
generally  accepted  theories,  regimes,  or  even  terminology.  Thus  even 
comparisons  of  results  and  statements  of  various  authors  are 
conplicatad.  But  when  conditions  are  properly  identified  and 
nomenclature  difficulties  are  overcome,  the  framework  to  be  presented  in 


this  work  is  found  to  be  coopatible  with  much  of  the  published  work. 


1.1  The  Jet  Breakup  Regimes 

Zf  ell  other  parameters  are  kept  constant,  the  jet  velocity  becomes  a 
convenient  quantity  to  introduce  various  regimes. 

Grant  and  Middleman*  reviewed  the  behaviour  of  low  speed  jets  and 
reported  the  results  in  the  form  of  a  breakup  curve  (Pig.  1)  which 
describes  the  unbroken  length  of  the  jet,  Li,  as  a  function  of  the  jet 
velocity,  U.  Once  a  jet  is  forsMd  (point  C,  Pig.  1}  the  jet  breakup 
length  at  first  increases  linearly  with  increasing  jet  velocity. 
Thereafter  it  reaches  a  maximum  (point  E)  and  then  decreases.  These 
first  two  breakup  regimes  are  reasonably  well  understood,  as  will  be 
seen  below,  and  here  are  called  the  Rayleigh  (O))  and  first  wind-induced 
breakup  (EP)  regimes.  A  feature  of  breakup  in  these  two  regimes  is  that 
drops  are  pinched  off  from  the  end  the  jet  and  their  sizes  are 
comparable  to  that  of  the  jet  (see  also  Pigs.  5a  and  5b}. 

Por  higher  velocity  jets  beyond  the  point  P  there  still  rmaains  some 
confusion  over  the  true  shape  of  the  breakup  curve.  Haenlein*  reported 
that  the  breakup  length  remains  constant,  or  decreases  slightly,  with 
increasing  velocity  (curve  PG,  Pig.  1)  and  then  it  abruptly  reduces  to 
near  zero  beyond  point  G.  This  suggests  the  existence  of  at  least  two 
more  breakup  regimes,  each  causing  new  features  in  the  breakup  curve. 
However,  McCarthy  and  Malloy*  and  Grant  and  Middleman*  report  that  the 
breakup  length  initially  increases  (curve  PH). 

It  should  be  noted  that  the  definition  and  measurement  of  the  intact 
length  becomes  increasingly  difficult  as  U  is  increased,  as  pointed  out 
also  by  Grant  and  Middleman*.  At  sufficiently  high  velocities,  the  jet 
surface  is  disrupted  prior  to  the  breakup  of  the  jet  core  and  the  use  of 
only  one  breakup  length  is  no  longer  a  cooqplete  measure  of  the  jet 


itabllity.  Thus  we  distinguish  between  the  intact'Surface  length, 
and  the  intact-core  length,  L,.  Zn  the  Kayleigh  and  first  wind-induced 
regimes  the  jet  breaks  up  simultaneously  over  the  entire  cross-section 
and  the  two  lengths  coincide.  Zn  the  second  wind-induced  and  atomization 
regimes  which  will  be  described  next,  the  disruption  starts  at  the  jet 
surface  and  eventually  reaches  the  jet  axis  so  that  at  least  two  lengths 
are  necessary  to  identify  the  gross  features  of  the  breakup. 

We  call  the  atomization  regime  that  regime  in  which  the  intact-surface 
length  is  zero  (but  the  intact-core  length  is  not  necessarily  zero)  so 
that  Lj  versus  U  follows  the  trend  suggested  by  Haenlein*,  viz.  the 
breakup  length  goes  to  zero  beyond  the  point  G  (H)  in  Pig.  1.  Here  the 
jet  surface  appears  to  break  up  immediately  at  the  nozzle  exit  and  drops 
are  formed  that  are  much  smaller  than  the  nozzle  diameter  (see  Pig.  5d). 
Although  we  find  this  definition  of  the  atomization  regime  unequivocal 
and  useful,  we  must  point  out  that  the  term  atomization  has  also  been 
used  by  other  authors  in  a  variety  of  different  contexts.  This  regime 
is  of  interest  in  many  fuel  injection  applications. 

Por  breakup  in  the  second  wind-induced  regime  (curve  PG  or  PH,  Pig.  1), 
both  the  intact-surface  length  and  the  intact-core  length  are  finite  and 
drops  are  formed  with  sizes  also  much  smaller  than  the  nozzle  diameter 
(see  Pig.  5c) . 

To  introduce  the  four  regimes  we  have  used  the  injection  velocity  as  a 
parameter  and  have  kept  all  other  parameters  constant.  Actually 
dimensionless  numbers  separate  the  various  regimes  as  will  be  discussed 
in  the  following  sections.  Indeed  the  high  velocity  jets  in  jet  cutting 
applications,  for  example,  exhibit  long  intact-surface  lengths  even 
though  they  are  faster  than  atomizing  fuel  jets,  but  the  two  families  of 
jets  differ  in  nozzle  design,  gas  density,  surface  tension  and  liquid 
viscosity. 

Attempts  have  been  made  by  various  authors  to  offer  criteria  with  which 
to  demarcate  breakup  regimes.  Por  example,  Kiesse^  correlated  breakup 


regime  data  and  prasentad  tha  rasults  in  a  form  suggastad  by  Ohnasorge* 
as  shown  In  Fig.  2.  Tha  boundarias  of  tha  ragimas  ara  raprasantad  by 
oblique  -straight  lines  on  a  graph  of  In  Z  (a  function  of  the  physical 
properties  of  the  liquid  and  tha  nozzle  diaoater  alone)  versus  In  Re, . 
Unfortunately,  this  method  of  correlation  does  not  include  the  effect  of 
the  initial  state  of  the  jet  (that  is  influenced,  for  axanple,  by  the 
nozzle  design)  nor  the  effect  of  tha  ambient  gas  density  (pressure), 
which  according  to  Torda*  modifies  tha  graph  as  shown  by  tha  dashed 
lines.  These  modified  boundary  curves  show,  for  azaa^la,  that 
atomization  can  be  achieved  at  lower  injection  velocities  by  injecting 
into  a  c«npressed  gas. 

The  effect  of  ambient  gas  density  on  jet  breakup  regimes  was  discussed 
by  Ranz**  who  argued  that  the  Weber  number  We,  should  be  a  controlling 
parameter.  He  offered  the  criterion  He,  >  13  for  the  onset  of 
atomization.  However,  it  should  be  noted  again  that  his  definition  of 
the  term  atomization  differs  from  ours  ~  he  does  not  refer  to  the  state 
of  the  jet  itself  but  Instead  refers  to  the  process  of  disintegration  of 
already  formed  droplets  during  their  flight  within  a  spray.  He  argued 
that  the  criteria  for  the  formation  and  the  subsequent  further  breakup 
of  the  droplets  should  be  the  same  since  when  the  inertial  stresses 
developed  by  the  surrounding  gas  exceed  the  surface  tension  stresses 
opposing  the  deformation  sufficiently,  the  liquid  drop  (or  ligament  in 
the  formation  process)  will  subdivide  into  smaller  units.  But  a  Weber 
number  correlation  by  itself  is  still  incomplete  since  now  the  liquid 
viscosity  is  not  accounted  for. 

This  latter  objection  can  be  removed  (at  least  conceptually)  by 
combining  Pig.  2  with  a  gas  density  parameter  as  is  shown  qualitatively 
in  Fig.  3,  taken  from  Reitz**.  Still,  different  sets  of  surfaces  should 
be  given  in  Pig.  3  to  account  for  the  effects  of  different  initial 
states  of  the  jet.  Thus  coopletely  satisfactory  correlations  for  the 
regime  boundaries  are  not  yet  available.  Indeed  many  authors  do  not 
even  distinguish  between  the  two  wind-induced  regimes  as  can  be  seen  in 
Fig.  2.  Finally,  as  previously  stated,  we  feel  that  part  of  the 
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difficulty  in  interpreting  published  results  is  due  to  a  lack  of  agreed 
upon  terminology  in  the  field,  and  insufficient  characterization  of  the 
injection  system  and  the  nozzle  geometry. 
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2.  LIHEAR  STABILITY  AKALYSIS 


W«  considsr  a  cylindrical  liquid  jat  issuing  from  a  circular  orifice 
into  a  stationary  incompressible  gas.  7he  stability  of  the  liquid 
surface  to  perturbations  is  axaminad  using  a  first  order  linear  theory 
which  ultimately  leads  to  a  dispersion  aquation,  Eq.  (5)  below.  This 
aquation  relates  the  growth  rate,  u,  of  an  initial  perturbation  of 
infinitesimal  anqplitude,  i),,  to  its  wavelength  X  (wavenumber  k«2ir/X}. 
The  relationship  also  includes  the  physical  and  dynamical  parameters  of 
the  liquid  jet  and  the  surrounding  gas.  The  present  treatment  will  be 
seen  (in  Section  3)  to  unify  the  results  of  Levich*>,  Stirling  and 
Sleicher^*  and  other  authors  who  have  treated  individual  jet  breakup 
regimes . 


The  column  of  liquid  is  assumed  to  be  infinite  in  the  axial  direction 
and  a  cylindrical  polar  coordinate  system  is  used  which  moves  with  the 
jet  velocity,  U.  Zn^sed  on  the  Initially  steady  motion  is  an 
infinitesimal  axisymmetric  surface  displacement,  one  Fourier  component 
of  which,  has  the  form 


The  linearized  hydrodynamical  equations  are 
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where  u^,  v^  and  Pi  are  small  axisymmetric  velocities  and  pressure. 
With  the  assumption  that  t7«a,  the  kinematic,  tangential  stress  and 
normal  stress  equations  are  to  first  order 


(2) 
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The  inertial  effects  of  the  gas  enter  through  the  gas  pressure  p, .  T1 
is  found  from  the  linearized  inviscid  equations  of  inotion  for  the  gas 

1^  ,  .  i_|E. 
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Pi3r 


where  the  mean  gas  motion  above  the  liquid  surface  is  given  by  U(r). 

The  boundary  conditions  are 

•  It  “  If  *t  r  a  a 

Ua»  Va,  Pa  *♦  0  aS  r  •  . 

Equations  (1)  are  solved  by  introducing  a  stream  function  ^ar  and  a 
velocity  potential  and  by  seeking  wave  solutions  of  the  form 

♦a  -  ♦a<r)  and  ^^a  •  *j(r) 

Solutions  free  from  singularities  on  the  axis  r^O  are  found  to  be  4>a 
Ca  Xa(kr)  and  <l>a  ■  Ca  rlaC-^r),  where  Ca  and  Cj  are  arbitrary  constants, 
and  the  liquid  pressure  can  then  be  found  from  the  relation 

Pa““Pi0^»/0t. 

For  the  gas  flow,  Eqs.  (3)  can  be  simplified  by  defining  a  stream 
function  ^(■(U-iw/k)n£(r} .  This  leads  to  an  Orr-Sossnerfeld  equation 

Hi  ♦  .  1)  511Z£)  .  k>f  -  0 

dr*  ^U-Wk  dr  ^  ^  ^ 

with  f(rBa)  ■  1,  and  f(r-*o}  ■  0.  The  equation  for  the  gas  pressure  is 
Pa  -  -P,»?(U-i|  )*  (fp  -  “)  •  < 


K  -  >  ~  -  ■  •  -  •  > 


H«re  tht  arbitrary  constant  ot  Integration  has  been  set  equal  to  zero. 
X£  the  gas  velocity  profile  U(r)  is  known,  the  gas  pressure  at  the  jet 
surface  can  be  determined  from  Eq.  (4)  for  use  in  Eq.  (2).  For  the 
special  case  of  slip  at  the  gas-liquid  interface,  U(r}»UBconstant,  and 
the  gas  surface  pressure  is 

...  K,(ka) 

p.  • -p.eu-ij)*!!,  . 


Finally,  substituting  these  relationships  into  Eqs.  (2)  yields 


(5a) 


which  is  the  governing  dispersion  relationship.  Equation  (5a)  may  for 
brevity  be  written  in  nondimensional  form  as 

fi*  *  2  Z  k*a*F»  0  -  kad  -  k*a*)  F*  ♦  We,  k*a*F,  (5b) 

where  0’»w7^S*7a,  2»M,//pIod  and  We,*p,U»d/o  and  the  P's  are 
dimensionless  ratios  of  Bessel  functions  and  wave  numbers. 
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3.  BREAKUP  REGIMES  AND  THE  DISPERSIOK  EQUATION 


In  this  Section  we  attempt  to  relate  the  mechanisms  of  breakup  in  the 
four  jet  breakup  regimes  to  limiting  cases  of  the  linear  stability 
analysis  of  Section  2.  This  theory  offers  a  unified  approach  for  the 
organization  of  the  jet  breakup  phenomena  but  is  not  complete  mostly 
because  it  does  not  account  explicitly  for  different  initial  states  of 
the  jet. 


3.1  Rayleigh  Breakup  Regime 

Rayleigh*  *  made  substantial  contributions  to  the  understanding  of  the 
stability  of  low  speed  jets.  He  obtained  a  dispersion  equation  for  the 
growth  of  axisymmetric  surface  disturbances  by  equating  the  potential 
and  kinetic  energies  of  an  inviscid  jet.  With  the  hypothesis  that  that 
disturbance  with  the  maximum  growth  rate  would  lead  to  the  destruction 
of  the  jet,  he  also  obtained  an  expression  for  the  droplet  size, 
assuming  that  it  would  be  of  the  order  of  the  wavelength  of  this 
disturbance. 


For  the  special  case  2*0  and  We,  «  o  (Inviscid  liquid  jet  at  low 
velocity)  the  dispersion  Equation  (5a)  (cf.  also  Eq.  [5b])  becomes 
.  ,  _  e  k 


P»a» 


(1  -  k»  a*) 


Ii {ka} 
I,{kai 


(6) 


which  is  the  same  as  Rayleigh's  result.  This  equation  predicts  that  the 
jet  surface  is  unstable  for  all  wavenumbers  with  ka<l  and  the 
corresponding  wavegrowth  curve  is  given  in  Fig.  4. 


This  wavegrowth  curve  can  be  found  experimentally  by  vibrating  low  speed 
jets  at  various  frequencies  and  by  measuring  the  growth  rate  of  the 
axisynnetric  surface  oscillations.  The  corresponding  measurements  of 
Donnely  and  Glaberson**  given  in  Fig.  4  show  excellent  agreement  with 
the  first  order  theory.  Differentiating  Eq.  (6)  shows  that  the  maximum 
growth  rate  is 


mnd  if  the  initial  disturbance  i)«  of  the  most  unstable  wave  grows 
exponentially  to  a  magnitude  ‘a*  in  time  T,  it  follows  that  the  breakup 
length  of  the  jet  (the  position  of  the  point  of  droplet  formation)  then 
on  average  will  be 

Li  -  UT  «  U  In  (a/ij.)  /  (7) 

lu 


This  linear  dependence  of  on  U  was  seen  in  Fig.  1  for  low  velocity 
jets  in  the  Rayleigh  breakup  regime. 


The  parameter  In  (a/r),)  has  been  determined  experimentally.  It  lies  in 
the  range  11-16  **  but  it  has  been  fouxvl  to  be  weakly  related  to  the 
Ohnesorge  number  Z  The  theoretical  influence  of  the  liquid 
viscosity  is  found  by  retaining  the  term  involving  Z  in  Eq.  (5b).  For 


large  Z  (high  liquid  viscosity)  the  maximum  wavegrowth  rate  is 

.  A.  .  q  i  ir  «  1 

"m  2/2  V.a*^  1  ♦  3Z  *  a/(Z+2) 


at  k 


where  ka  has  been  assumed  to  be  small  and  the  Bessel  functions  and  their 
arguments  have  been  replaced  by  their  asymptotic  values.  This 
relationship  was  first  obtained  by  Weber**.  His  analysis  showed  that  the 
effect  of  the  increased  liquid  viscosity  is  to  move  the  most  unstable 
wave  to  longer  wavelengths  without  altering  the  value  of  the  stability 
boundary,  ka«l.  The  jet  breakup  agency  remains  the  destabilizing 
combination  of  surface  tension  and  Inertia  forces  on  the  jet. 


The  above  analysis  predicts  that  the  jet  breakup  yields  droplets  many 
nozzle  diameters  downstream  of  the  nozzle.  The  drop  diameters  are  larger 
than  that  of  the  jet  and  a  photograph  typical  of  jet  breakup  in  this 
regime  is  shown  in  Fig.  5a,  taken  from  Lee  and  Spencer**. 


To  estimate  the  droplet  size,  Rayleigh  assumed  that  all  of  the  liquid 
enclosed  within  the  wave  forms  the  volume  of  the  newly  created  drop  when 


the  surface  wave  amplitude  equals  the  jet  radius.  However,  as  pointed 
out  by  Nang**  when  the  wave  an^litude  becomes  comparable  to  the  jet 
radius,  the  surface  deformation  is  observed  to  be  non-sinusoidal  due  to 
nonlinear  effects.  Also,  mass  is  only  conserved  to  first  order  in  the 
first  order  stability  analysis.  Yuen**  and  Kayfey**  retained  higher 
order  terms  in  their  jet  stability  analysis  and  Rutland  and  Jameson** 
demonstrated  that  this  improved  theory  also  predicts  the  existence  of 
satellite  droplets  formed  between  the  primary  drops.  Interestingly, 
LaPrance**  showed  that  the  drop  sizes  are  not  influenced  by  the 
magnitude  of  the  initial  disturbances  to  the  jet. 


3.2  First  Wind-Induced  Breakup  Regime 


The  second  term  on  the  right  hand  side  of  Eqs.  (5)  becomes  important 
when  the  jet  velocity  (for  example)  is  increased,  in  this  case  the 
inertial  effects  of  the  surrounding  gas  can  no  longer  be  neglected  and 
the  Weber  number  Wej  becomes  a  controlling  parameter  in  the  dispersion 
equation.  Weber**  showed  that  the  effect  of  the  environment  on  the  jet 
is  to  enhance  the  growth  rate  of  disturbances,  leading  to  earlier 
breakup  of  the  jet.  He  obtained  the  result 

«*  ♦  3K,k*u  -  2^,(1  '  h*a*)  k*a* 

U*k*a»  K.(ka) 

Pi  2a*  Ki(ka) 


which  can  also  be  found  from  Eqs.  (5)  in  the  limit  ka  <  1. 

Computations  of  the  maximum  wave  growth  rate  now  show  that  the  jet 
breakup  length  equation,  Eq.  (7),  becomes  a  nonlinear  formula  which  must 
be  solved  numerically.  It  indeed  predicts  that  the  breakup  length  Lx 
decreases  with  Increasing  jet  velocity  (Fig.  1,  curve  EF).  However  the 
predicted  value  of  the  maximum  in  the  breakup  curve  (Fig.  1,  point  E) 
fails  to  agree  with  experimental  results,  since  Weber's  theory  is  found 
to  overestimate  the  aerodynamic  effect  of  the  gas.  This  has  led  many 
investigators  to  attempt  modifications  to  Eq.  (8). 


Fenn  and  Middleman**  argued  that  the  viscosity  of  the  gas  should  also  be 


considcrad  in  a  oora  cooplata  analysis.  The  affect  of  the  gas  viscosity 
enters  through  the  normal  and  tangential  stress  boundary  conditions  Eqs. 
(2)  and  was  neglected  by  Weber  (also  in  the  linear  stability  analysis  of 
Section  2).  The  inviscid  gas  approximation  implies  slip  at  the  gas 
liquid  interface. 

Benjamin’*  showed  that  if  the  gas  boundary  layer  is  thin  compared  to  the 
surface  wave  wavelength  (high  gas  Reynolds  numbers)  then  the  energy 
transmitted  between  the  gas  and  the  liquid  by  normal  stresses 
(aerodynamic  effect)  is  large  compared  to  that  transmitted  by  shear 
stresses  (viscous  effect)  ,  and  their  ratio  is  independent  of  the 
Reynolds  number.  However,  for  finite  Reynolds  numbers  the  magnitude  of 
the  fluctuating  pressure 'component  pj  which  is  in  phase  with  the  wave 
elevation  (the  part  responsible  for  energy  transfer)  becomes  reduced. 
Based  on  these  results,  Sterling  and  Sleicher”  introduced  an 
attenuation  coefficient  C,  smltiplying  p,  with  C.«  l-h(ka,Re2  ).  The 
function,  h,  increases  as  k  increases  and  decreases  to  zero  as  the 
Reynolds  number  becomes  large.  They  found  that  C  could  be  replaced  by  a 
constant,  noticing  that  the  wavenumber  and  the  Reynolds  number  move  in 
opposite  directions  as  the  aerodynamic  effects  increase. 

This  modification  of  Weber's  theory  with  C  «  0.175  multipying  the  second 
term  on  the  right  hand  side  of  Eqs.  (5)  agrees  well  with  experimental 
results.  The  numerical  results  of  Sterling  and  Sleicher**  then  show 
that  the  point  E  in  Pig.  1  is  reached  when  We,  ■  1.2  3.41Z***. 

In  this  regime  the  jet  breakup  still  occurs  many  nozzle  diameters 
downstream  of  the  nozzle  and  produces  droplets  whose  diameters  are  still 
comparable  to  that  of  the  jet  as  was  seen  in  Pig.  5a  for  jets  in  the 
Rayleigh  breakup  regime  (see  Pig.  5b).  The  breakup  is  still  due  to  the 
destabilizing  influence  of  surface  tension,  but  it  is  now  augmented  by 
the  aerodynamic  interaction  between  the  liquid  and  gas. 

Sterling  and  Sleicher”  also  pointed  out  that  relaxation  of  the  jet's 
exit  velocity  profile  to  a  uniform  flow  beyond  the  nozzle  can  influence 
the  jet  ;.xeakup  process.  In  their  analysis  (which  was  quoted 
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abova)  they  considarad  tha  case  whera  l^/l  Is  small  (1^  is  tha  jet 
proflla  ralaxatlon  langth).  For  large  1^/L  they  suggest  that  jet 
instability  could  be  enhanced  by  velocity  profile  rearrangement  effects 
but  they  offered  no  details  of  the  mechanism  by  which  this  would  occur. 
The  profile  relaxation  phenomenon  was  also  alluded  to  by  Grant  and 
Middleman*  who  found  a  dependence  of  Li  on  the  nozzle  passage  length.  In 
particular,  nozzles  of  short  length  produced  more  stable  jets  than  those 
produced  from  long  tubes.  They  argued  that  this  implies  a  coupling 
betveen  the  velocity  profile  and  the  mechanism  of  instability  in  this 
breakup  regime.  Parenthetically,  their  experiments  also  showed  that  the 
breakup  length  in  the  Rayleigh  regime  (Section  3.1)  is  not  influenced  by 
nozzle  design  details. 

Another  example  of  the  influence  of  the  initial  state  of  the  jet  is 
discussed  by  Phinney*'  who  proposed  that  liquid  turbulence  also  enhances 
the  jet  breakup  process.  He  reasoned  that  the  effect  of  the  jet 
turbulence  is  to  increase  the  initial  disturbance  level  n**  He  noted 
that,  even  in  the  abscence  of  aerodynamic  effects,  the  jet  breakup 
length  is  reduced  once  a  critical  value  of  the  jet  Reynolds  number  is 

reached.  Furthermore,  this  critical  Reynolds  number  is  of  the  same 
order  as  that  for  transition  to  turbulence  in  the  nozzle.  However  the 
influence  of  jet  turbulence  on  the  magnitude  of  the  initial  disturbance 
level  is  still  quantitatively  unclear. 


3.3  Second  Wind-Induced  Breakup  Regime 


With  further  increases  in  Wej  ,  Eq.  (8)  predicts  that  the  maximum  wave 
growth  rate  occurs  at  progressively  larger  wavenumbers  (shorter 
wavelengths).  An  inspection  of  Eq.  (8)  (or  Eqs.  [5]),  shows  that  the 
first  term  on  the  right  hand  side  changes  sign  at  ka*!,  after  which  the 
surface  tension  forces  oppose  the  breakup  process.  Jet  breakup  is  now 
due  to  the  unstable  growth  of  short  wavelength  surface  waves  (ka>l) 
which  are  induced  by  the  relative  motion  between  the  jet  and  the  ambient 
gas.  An  analysis  of  Eq.  (8)  shows  that  the  maximum  wave  growth  rate 
occurs  at  ka>l  when  We,aI2  for  inviscid  jets.  This  estimate  was  made 


using  the  numerical  results  of  Sterling  and  Sleicher**.  The  estimate 
also  agrees  well  with  the  experimentally  obtained  criterion  of  Ranz^*, 
Wea«13,  for  the  onset  of  short  wavelength  waves. 

An  expression  for  the  growth  rate  of  short  wavelength  surface  waves  was 
presented  by  Levich^*  and  Levich  and  Krylov* who  examined  Eqs.  (5)  in 
the  limit  ka  »  1  and  deduced,  neglecting  the  liquid  viscosity,  that 

-  (Pak*U*  -  ok*)/pa  .  (9) 


This  result  in^lies  the  existence  of  unstable  surface  waves  when  k  < 
p,U*/e.  The  maximum  groirth  rate  is  given  by 
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Equation  (9)  shows  that  the  dispersion  relation  Eq.  (5)  becomes 
independent  of  the  jet  radius  in  this  limit.  Consequently,  for  ka  »  1 
jet  curvature  effects  are  unimportant.  Similarly,  the  Weber  number  We, 
can  no  longer  appear  as  a  controlling  parameter. 

The  influence  of  the  liquid  viscosity  is  seen  by  retaining  the  second 
term  on  the  left  hand  side  of  Eqs.  (5).  In  the  limit  ka  this  reduces 
to 

(«  +  Ivjk*)*  ♦  ok*/p,  -  4i',*k*/(k*  ♦  u/vi) 

♦  (w  ♦  iUk)*p,/pj  ■  0  . 


This  result  is  identical  to  that  of  Taylor**  who  performed  an  analysis 
of  the  unstable  growth  of  2-dimensional  planar  surface  waves  due  to  the 
relative  motion  between  a  liquid  and  a  gas.  He  considered  the  limit  ka 
»  1  and,  assuming  p,  «  P|  he  found  that  the  wave  growth  rate  is 

w/kU  -  2  (  ^  K*xg(r,x)  . 

Pj 


The  function  g  is  a  correction  to  the  result  of  Levich**  which  now 
accounts  for  the  effect  of  the  liquid  viscosity.  It  is  shown  as  a 
function  of  the  new  parameter 


•nd  the  nondinensional  wavelength  x  >  from  Taylor's  work  in  Pig. 

6.  The  figure  shows  that  the  disturbance  growth  rate  increases  with 
increasing  T,  and  that  the  maximum  growth  rate  occurs  at  larger 
wavenumbers  (shorter  wavelengths)  as  T  increases. 


Taylor**  also  estimated  the  intact-core  length,  L,,  by  coeluting  the 
rate  at  which  droplets  remove  mass  from  the  liquid  core.  Here  the 
dr^let  sizes  were  assumed  to  be  proportional  to  the  unstable  surface 
wave  wavelengths.  This  analysis  gives  (see  also  Reitz  and  Bracco**) 

L,  -  B*  a  (  )V’f(r)  (10a) 

Pi 


where  Bi  is  a  constant  of  order  unity.  The  function  f(r)  corresponds  to 
the  maximum  wavegrowth  rates  of  Pig.  6,  and  it  is  shown  in  Pig.  7.  The 
intact-surface  length  Lj  can  be  estimated  using  similar  arguments  to 
those  which  led  to  the  development  of  Eq.  (7)  (where  the  jet  radius  'a* 
is  replaced  by  some  characteristic  wave  height  at  breakup,  say).  In 
this  case  we  find  that 

Li  -  B,  L,  /  We,  (10b) 


where  B,  is  another  constant  which  would  be  dependent  weakly 
(logarithmically)  on  drop  size. 


Equation  (10a)  predicts  that  the  intact-core  breakup  length  L,  remains 
constant  with  increasing  jet  velocity  until  the  parameter  r  becomes 
small  (i.e.  for  high  velocity  jets).  L,  then  decreases  with  increasing 
jet  velocity.  The  intact-surface  length  L,  from  Eq.  (10b)  is  predicted 
to  decrease  as  the  inverse  square  of  the  velocity  and  then  to  decrease 
even  faster  when  T  becomes  small.  This  prediction  is  not  inconsistent 
with  the  form  of  breakup  curve  of  Haelein*  for  jets  in  the  second 
wind-induced  regime  and,  possibly,  even  in  the  atomization  regime  once 
L,  becomes  of  the  order  of  the  wavelength  of  the  surface  waves. 
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A  photograph  typical  of  jat  braakup  in  this  raglaie  Is  shown  in  Pig.  5c 
fron  Raitz**.  Tha  photograph  shows  that  tha  jat  braakup  starts  some 
dlstanca  downstraan  of  tha  nozzle  exit  and  yields  droplets  whose  average 
diameters  are  much  lass  than  tha  jet  diameter.  Droplet  formation  results 
from  tha  unstable  growth  of  short  wavelength  surface  waves  on  the  jet. 
This  wave  growth  is  aused  by  the  relative  motion  between  the  jet  and 
the  ambient  gas,  and  surface  tension  forces  oppose  the  wave  growth 
process. 


3.4  Atomization  Regime 

In  this  regime  the  breakup  appears  to  comaence  at  the  nozzle  exit  (see 
Pig.  5d).  The  spray  takes  the  form  of  a  cone  with  its  vertex  within  the 
nozzle.  Various  authors  have  suggested  possible  jet  breakup  agencies 
for  jets  in  the  atonization  regime  and  some  of  these  are  considered 
below,  but  a  conplete  and  tested  theory  is  not  yet  available. 

Experiments  were  made  to  study  atcznization  in  Reitz  and  Bracco*  ^ *  and 
Wu  et  al . ^ .  The  range  of  conditions  of  these  studies  has  been 
extensive  and  includes  the  operating  conditions  of  fuel  injection 
systems  in  Diesel  and  stratified  charge  internal  combustion  engines. 

The  experiments  were  made  under  steady  conditions  with  injections  into  a 
semi'infinite  gas.  The  test  conditions  include:  Constant  liquid 
Injection  pressures  in  the  ranges  3.4-17.0  MPa*^  and  10.8-90.5  MPa**; 
Constant  gas  pressures  in  the  range  0. 1-4.1  MPa  with  air,  nitrogen, 
helium  and  xenon  (different  molecular  weights  to  isolate  effects  of  gas 
density  and  pressure**};  Water  and  water  glycerol  injections  (10* 
range  in  liquid  viscosity**}.  Hexane,  water  and  tetradecane  (factor  of 
10  in  viscosity,  4  in  surface  tension,  1.5  in  liquid  density**.  Pentane, 
hexane  and  ethanol  (factor  of  3  surface  tension,  1.5  liquid  density*;  21 
nozzles:  sharp  edge  inlet  tube  nozzles  (length  to  diameter  ratio  range 
0.5-65.0  (diameter  0.35  mm},  rounded  inlet**  and  rounded  exit  nozzles 
and  cavitation  free  nozzles,  and  a  factor  of  2.5  in  nozzle  exit 
diameter**’*;  Liquid  tenperature:  room  temperature**'**:  100-200  C*. 


certain  trands  wara  found  froo  photographs  of  jat  braakup  In  the 
axparlDants.  For  aaanple,  the  spray  angle  (divergence  angle  of  the  jet) 
was  found  to  increase  with  increasing  (isothermal)  gas  compression. 
Moreover,  it  was  established  that  this  is  due  to  increases  in  the  gas 
density,  not  pressure**.  The  spray  angle  was  also  found  to  increase  with 
decreasing  liquid  viscosity. 

Another  trend  concerns  the  boundary  between  the  atomization  and  the 
second  wind-induced  regimes.  It  was  found  that  the  breakup  starts 
progressively  closer  to  the  nozzle  exit  as  the  gas  density  is  increased, 
until  it  reaches  the  exit  with  no  evidence  of  an  abrupt  change.  This 
trend  is  also  shown  in  Pig.  8  in  which  the  measured  spray  angle  is 
plotted  against  the  gas-liquid  density  ratio  for  sharp  edge  inlet 
nozzles  with  a  length  to  diameter  ratio  of  4.0.  The  solid  data  points 
are  in  the  second  wind-induced  regime  -  jets  intact  before  diverging. 

The  open  points  show  atomizing  jets  -  divergence  begins  at  the  nozzle 
exit.  It  was  also  found  that  atomization  is  reached  once  the  liquid 
viscosity  is  decreased  below  a  certain  level,  again  with  no  abrupt 
change  in  the  appearance  of  the  jet. 

Other  results  were:  The  spray  angle  decreases  with  increased  nozzle 
passage  length  for  nozzle  length  to  diameters  greater  than  10  or  20. 

(For  shorter  nozzles  there  is  more  scatter  in  the  results  and  the  trends 
are  not  yet  fully  established,  Hu  et  al.**);  For  the  same  length, 
roundi^  inlet  nozzles  produce  less  divergent  jets  than  sharp  edge  inlet 
nozzles;  Atomization  commences  at  different  gas  density  and  liquid 
viscosity  levels  as  the  nozzle  design  is  changed. 

With  these  results  it  is  possible  to  examine  previously  proposed 
theories  for  atomization  in  detail.  For  exanf>le,  DeJuhasz**  and 
Schwletzer**  proposed  that  liquid  turbulence  causes  atomization.  But  if 
pipe  turbulence  were  the  only  mechanism,  turbulent  jets  (from  the 
nozzles  with  large  length  to  diameter  ratios**)  would  have  been  the  most 
unstable  flows  -  contrary  to  the  experiments.  Similarly,  cavitation 
phenomena  were  proposed  Bergwerk**  to  lead  to  atomization.  But  jets  were 


found  to  atonize  ovon  when  the  cavitation  free  nozzles  of  Wu  et  al.** 
were  used. 

Zn  fact,  an  evaluation  of  other  proposed  atomization  mechanisms  has 
revealed  that  none  of  the  theories  taken  alone,  is  able  to  explain  the 
results  fully**.  These  theories  include  proposals  that  atomization  is 
caused  by:  Aerodynamic  surface  wave  growth**'**  -  the  results  would  be 
Independent  of  the  nozzle  geometry;  Rearrangement  of  the  jet's  velocity 
profile**  >  the  high  viscosity  jets  would  be  the  most  unstable;  Liquid 
supply  pressure  oscillations**  -  atomization  would  not  have  occurred 
since  the  pressure  was  constant  in  the  experiments;  and  Wall  boundary 
layer  velocity  profile  relaxation**  **  atomization  would  have  been 
Independent  of  the  gas  density;  all  of  these  contrary  to  the 
experiments. 

However,  the  aerodynamic  surface  wave  growth  theory  was  found  to  predict 
many  of  the  trends  in  tests  with  a  given  nozzle.  Thus  it  is  useful  to 
consider  its  predictions  in  more  detail.  In  this  case,  the  appropriate 
limit  of  the  dispersion  Equation  (5)  is  ka  **  •,  as  in  the  second 
wind* induced  regime  above.  Ranz**  argued  that  the  spreading  angle  of 
the  atomizing  jet  could  be  predicted  by  combining  the  radial  velocity  of 
the  fastest  growing  of  the  unstable  surface  waves  with  the  axial 
injection  velocity: 

tan  («/2)  -  5  -  -  4x  (pj/p^Kf  (D/A  (11) 

where  the  proportioiMlity  constant  A  is  obtained  from  experiment. 

In  Eq.  (11)  the  spray  angle  increases  with  increasing  gas  density  and 
decreases  with  increasing  liquid  viscosity  and  increasing  velocity. 

This  is  consistent  with  tests  with  a  given  nozzle.  Notice  also  from 
Pig.  7  that  for  T  >  1  the  function  f(r}  becomes  asymptotically  equal  to 
✓3/6.  Equation  (11)  then  predicts  that  the  spray  angle  depends  only  on 
the  gas-liquid  density  ratio  -  which  is  surprising  considering  the  many 
parameters  that  could  effect  it.  This  behaviour  is  also  generally  borne 
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out  by  the  Beesurenentt  with  e  given  nozzle 


Renz«*  determined  that  A  »  16  or  20,  but  he  pointed  out  that  the  data  of 
Schweitzer* «  give  A  »  3.  Results  given  in  Reitz  and  Bracco**  indicate 
that  the  predicted  variation  of  the  spray  angle  with  gas  density  is 
followed  if  A  has  a  different  value  for  each  different  nozzle.  Also 
other  results  show  that  agreement  is  found  with  respect  to  liquid 
viscosity  variations  by  using  (for  each  nozzle)  the  same  value  of  A  as 
that  obtained  from  the  gas  density  best  fit.  Figure  9  shows  spray  angle 
measurements  and  the  predictions  of  Eq.  (11)  with  respect  to  gas 
density  and  jet  velocity  variations.  Again  the  measurements  follow  the 
predicted  trend,  but  the  results  do  exhibit  a  mild  opposite  trend  at  low 
gas  densities,  indicating  that  the  theory  only  complies  with  the  most 
pronounced,  and  practically  important,  of  the  measured  trends. 

There  is  additional  evidence  in  support  of  the  aerodynamic  theory.  For 
exaiqile,  the  data  of  Hiroyasu  at  al.**  agree  with  the  prediction  of  the 
intact-core  length  of  Eq.  (10a)  with  respect  to  injection  velocity  and 
gas  density  changes.  Their  experiments  were  based  on  measurements  of 
electrical  resistance  between  the  nozzle  and  a  screen  that  could  be 
moved  axially  within  the  spray;  thus  detecting  any  continuous  liquid 
connection  between  the  nozzle  and  the  screen.  However,  a  connected 
ligament  could  give  the  same  signal  as  a  solid  liquid  core  free  of  gas. 
Consequently  there  is  still  some  uncertainty  about  the  structure  of  the 
core. 

Drop  sizes  have  been  measured  at  the  edge  of  the  spray  in  the  vicinity 
of  the  nozzle  exit**.  In  the  surface  wave  growth  theory  the  drop  size 
would  be  related  to  the  wavelength  of  the  unstable  surface  waves  i.e. 

r»  c  2wcx/p,U*  (12) 

where  c  is  a  constant  of  order  unity.  The  predictions  of  Eq.  (12)  with 
regard  to  liquid  properties,  jet  velocity  and  lack  of  sensitivity  to 
nozzle  design  are  found  to  be  in  agreement  with  experiment.  However  a 
discrepancy  has  been  found  with  the  deptfuSence  on  gas  density  but  it 


•ppaars  that  a  raasonable  axplanatlon  for  this  dlsagraament  exists*'. 


Based  on  the  definition  that  atooization  occurs  when  the  intact  surface 
length  approaches  zero,  Eq.  (10a)  can  be  used  to  predict  that 
atomization  occurs  if: 

(Pi/P,)'^*<  K  for  r  >  1 

and  (  )*'^*<  K  for  r  <  1  (13) 

Pa 

where  the  constant,  K,  depends  on  the  nozzle  geometry.  For  the  nozzles 
in  Pig.  8  (length  to  diameter  ratio  4),  K  is  found  empirically  to  be 
equal  to  9.2. 

The  stability  analysis  does  not  include  the  details  of  the  flow  field 
within  the  nozzle.  This  shortcoming  was  also  mentioned  for  the  first 
and  second  wind-induced  breakup  regimes  in  Sections  3.2  and  3.3.  Thus, 
it  appears  that  the  stability  analysis  can  explain  many  of  the 
experimental  results  in  the  atomization  regime  (and  also  in  the  second 
wind-induced  regime)  if  it  is  supplemented  by  additional  information 
pertaining  to  the  initial  state  of  the  jet  since  such  information  is  not 
included  explicitly  in  the  theory.  In  particular,  the  nozzle  geometry 
obviously  effects  the  initial  state  of  the  jet.  The  simplest  way  of 
accounting  for  the  initial  state  of  the  jet  is  through  the  magnitude  of 
the  initial  perturbation  ?}•  which  could  have  a  different  value  for 
different  nozzle  geometries,  for  example.  However,  even  if  this  is 
eventually  found  to  be  a  sufficient  modification,  the  magnitude  of  the 
perturbation  is  known  at  best  only  enpirically.  All  that  can  be  said  at 
present  is  that  ?;•  depeiuls  in  some  complex  manner  on  the  details  of  the 
initial  jot  flowfield  that  in  turn  is  influenced  by  the  nozzle  geometry 
and  possibly  by  other  parameters,  such  as  the  flowfield  just  upstream  of 
the  nozzle. 

However,  a  physical  picture  of  the  atomization  process  can  now  be 
proposed  which  is  consistent  with  the  available  data.  The  surface  of 
the  liquid  jet  emerges  from  the  nozzle  already  perturbed  by  events  that 
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occur  within  the  nozzl*  itself  and  ere  effected  by  its  geometry.  The 
perturbations  are  rapidly  and  selectively  amplified  by  aerodynamic 
interaction  with  the  gas  until  the  outer  surface  of  the  jet  breaks  into 
drops.  The  size  of  the  drops  and  the  intact-surface  length  is  much 
smaller  than  the  diameter  of  the  nozzle.  The  depth  from  the  surface  of 
the  jet  to  which  the  above  drop  formation  mechanism  would  apply  is  not 
known.  But  the  core  eventually  breaks  up  too  since  only  isolated  drops 
are  found  far  downstream. 

Other  aspects  of  atomization  still  remain  unresolved  besides  the 
influence  of  the  nozzle.  Not  predicted  by  the  stability  theory,  and 
therefore  unknown,  are  the  size  and  size  distribution  of  the  unstable 
waves  at  the  moment  of  breakup  and  also  the  time  between  successive 
ruptures.  Away  from  the  nozzle  exit,  as  the  generating  surface 
regresses  towards  the  axis  of  the  jet,  there  are  questions  as  to  what 
gas  velocities  are  seen  by  the  liquid  surface.  The  velocity  of  the 
entrained  gas  certainly  approaches  that  of  the  liquid  surface.  Thus  the 
breakup  process  should  be  coupled  with  the  two-phase  flow  field  that 
exists  between  the  presumed  intact-core  and  the  unperturbed  outer  gas. 

As  the  relative  velocity  between  the  liquid  and  gas  decreases  inside  the 
jet,  larger  drops  or  ligaments  or  blobs  should  be  formed;  just  as  larger 
drops  are  found  when  the  injection  velocity  is  decreased,  as  in  going 
from  Pig.  5d  to  Fig.  5a. 

An  additional  factor  to  be  considered  is  coalescence  of  the  liquid 
fragments  which  can  be  expected  where  locally  large  values  of  the  liquid 
volume  fraction  exist.  Xn  fact,  the  net  outcome  of  atomization  may  be 
the  result  of  a  small  difference  between  large  droplet  formation  and 
coalescence  rates**.  Thus  it  is  clear  that  Bqs.  (10)  to  (13)  may 
provide  some  information  about  the  outcome  of  the  breakup,  but  in  no  way 
do  they  give  all  that  is  necessary. 
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wave  number  2ir/X 
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jtt  bTMXup  time 
axial  velocity  component 

jet  exit  velocity  (averaged  over  jet  cross-section) 

radial  velocity  component 

Weber  number  pU*d/o 

dimensionless  wave  length  pjUVoh 

axial  coordinate 

Ohnesorge  number  Pi/^(Piod) 

dimensionless  wavegrowth  /o) 

Taylor  parameter  ■PiO*/(p»i»i*U*) 
effective  injection  pressure  Pj-Pa 
surface  wave  amplitude 
jet  divergence  or  spray  angle 
wave  length 
kinematic  viscosity 
dynamic  viscosity  ■  p/p 
fluid  density 

surface  tension  coefficient 
velocity  potential 
stream  function 
wave  growth  rate  "ftCu) 


Subscripts 

1,L  liquid  phase 
2>g  gas  phase 

m  maximum  value 


FIGURES 


1.  Jet  surface  breakup  length  Lj  as  a  function  of  jet  velocity:  ABC 
Drip  flow;  CD  Rayleigh  breakup  regime;  EF  First  wind-induced  breakup 
regime;  PG  (PH)  Second  wind-induced  breakup  regime;  Beyond  G  (H) 
Atomization  regime. 

2.  Jet  breakup  regime  boundaries  of  Miesse\  Ohnesorge*  and  Torda*. 

3.  Schematic  chart  of  influence  of  gas  density  on  breakup  regime 
boundaries «  Reitz  ^  * . 

4.  Wave  growth  rate  with  wave  number  for  jets  in  Rayleigh  breakup 
regime.  Line:  theory  Eq.  (6).  Symbols:  measured  wave  growth  rates  of 
Donnely  and  Glaberson^*. 

5.  Ezan^les  of  jets  in  the  four  breakup  regimes  a.)  Jet  breakup 

in  the  Rayleigh  regime.  Droplet  sizes  of  the  order  of  the  jet  diameter 
and  breakup  occurs  many  nozzle  diameters  downstream  of  the  nozzle,  b.) 
Jet  breakup  in  the  first  wind-induced  regime.  Droplet  sizes  are  still 
of  the  order  of  the  jet  diameter  and  breakup  occurs  many  nozzle 
diameters  downstream  of  the  nozzle,  c.)  Jet  breakup  in  the  second 
wind-induced  regime.  Droplet  sizes  much  smaller  than  the  jet  diameter 
and  the  breakup  starts  some  distance  downstream  of  the  nozzle,  d.)  Jet 
breakup  in  the  atomization  regime.  Droplet  sizes  much  smaller  than  the 
jet  diameter  and  the  breakup  starts  at  the  nozzle  exit. 

6.  Theoretical  wave  growth  rate  as  a  function  of  dimensionless  wave 
length  parameter,  p,U*/ok  for  jets  in  the  second  wind-induced  breakup 
regime,  Taylor*  * . 

7.  DepeixSence  of  the  maximum  growth  rates  in  Fig.  6  on  Taylor's 
parameter  r  ■  (pi/Pi){Rei/Wei)*. 

6.  Measured  spray  angle  versus  gas-liquid  density  ratio  for  nozzle 
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length  to  diameter  ratio  of  4  nozzle  diameter  0.34  m.  Nozzle 

dP>10. 0-14.0  MPa.  Liquid:  water.  Gas:  nitrogen.  Solid  symbols:  jets 
Intact  before  diverging  (second  wind-induced  regime).  Open  symbols:  jet 
divergence  starts  at  the  nozzle  exit  (atomization  regime).  Partially 
solid  symbols:  Marginal  breakup  at  the  nozzle  exit.  Line:  aerod^nnamic 
theory  prediction  Eq.  (11)  with  A  >  4.0. 

9.  Measured  spray  angle  versus  liquid  jet  velocity  for  different 
gas-llquld  density  ratios’*.  Liquid:  n-hexane.  Gas:  nitrogen.  Straight 
tube  nozzle  length  to  diameter  ratio  4.  Open  symbols:  data  at  AP  >  15.3, 
38.0,  64.9  and  91.8  MPa  (U  -  1.5,  2.3,  3.0  and  3.6xl0«  cm/s, 
respectively).  Closed  symbols:  repeated  measurements  at  AP«15.3  MPa  made 
after  the  high  Injection  pressure  (AP-91.8  MPa)  runs. 
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LDV  Measurements  of  Drop  Velocity  in  Diesel-Type  Sprays 

K.-J.  Wu,*  D.  A.  Santavicca.t  and  F.  V.  BraccoJ 
Princeton  University,  Princeton,  f^Jew  Jersey 
and 

A.  Coghe§ 

Centro  Ricerche  Propulsione  e  Energetica  CNR,  Milan,  Italy 

Axial  and  radial  componenls  of  ibc  drop  vrlodly  were  measured  by  laser  Doppler  velocimelry  (LDV)  wllbln 
n-hexane-lnlo<nllrogen  sprays  from  single-hole  cyliiMittcal  nozdes  al  room  lemperalurc.  The  gas-to-liquid 
density  ratio,  the  injection  velocity,  and  the  nozzle  knglh  and  diameter  were  varied.  The  LDV  system  Included 
an  argon-ion  laser,  dual  beam  LDV  optics  with  frequency  shifting  and  90-deg  collection,  a  counter  processor, 
and  minicomputer  data  acquisition.  The  total  erroi  in  the  measured  mean  nnd  fluctuation  drop  velocities  was 
less  than  lOVo  from  the  centerline  to  the  half-radius  (half-tbe-wldth  al  balf-lhe-deplh)  nnd  larger  and  more 
uncertain  beyond  It.  Velocity  bias  was  the  most  difficull  error  to  quantify.  It  Is  found  that  beyond  300  nozzle 
diameters  from  the  nozzle  so  much  ambient  gas  has  been  entrained  by  the  drops  that  the  subsequent  struclure  of 
the  spray  is  dominated  b>  the  entrained  ambient  gas,  and  the  fully  developed  Incompressible  Jet  structure  and 
drop-gas  equilibrium  are  being  approached.  This  conclusion  is  supported  by  all  the  measured  drop  velocity 
parameters:  Jet  half  radius,  centerline  velocity  decay,  axial  mean  velocily  distribution,  axial  and  radial  velocity 
fluctuation  distributions,  nnd  Independence  of  drop  velocity  on  drop  size.  Large  uncertainties  about  the 
magnitude  of  errors  at  the  outer  edges  of  Jets  using  both  laser  Doppler  velocimelry  and  hot-wire  nnemometry 
suggest  that  this  region  is  still  poorly  characterized  even  for  incompressible  Jets. 


Nomcnctftture 

C.C,,C,  =  constants  of  Eqs.  (I )  and  (2) 

d  =  nozzle  orifice  diameter,  yim 

=natness=  I',"  / \V;‘ ]> . 

L  =  nozzle  passage  length,  ytm 

p,  =  chamber  gas  pressure,  M Pa 

r  =  radial  coordinate,  cm 

r^s  =  half-radius  (half-the-width  at  half-the-depih)  of 

the  drop  axial  velocity  profile,  cm 
RFA,  ,  =  relative  fluctuation  amplitude  (standard 

deviation)  of  the  axial  velocily  of  the  drops, 

=  s'Vj/V, 

=ske«ness= 

U,,  U,  =  radial  and  axial  components  of  the  gas 
velocity,  m/s 

U„  ci  ~  velocity  of  the  gas  at  the  center- 
line  of  the  jet,  m/s 

l'„  F,  =  radial  and  axial  componenls  of  the 
drop  velocity,  m/s 
=  mass  mean  injection  velocity,  m/s 
P,  =  axial  velocity  of  the  liquid  at  the 

centerline  of  the  spray,  m/s 
X,Y,Z  =  spatial  coordinates  (X  is  from  the 

nozzle  exit,  see  Fig.  3),  cm 
Xg  =  virtual  origin,  cm 

Ap  -  effective  injection  pressure,  MPa 

6  =  gas  volume  fraction 

p,  =  liquid  dynamic  viscosity,  kg/m  -  s 

Pj,  =  chamber  gas  density,  kg/m’ 

p,  =  liquid  density,  kg/ m' 

o,  =  liquid  surface  tension  coefficient,  kg/s* 
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Superscripts 

(■)  =  mean  value 

(  )'  =  fluctuating  component 


Introduction 

HE  breakup  of  liquid  jets  is  achieved  through  a  large 
variety  of  atomizers  for  an  even  larger  variety  of 
applications. Its  purpose  is  to  increase  the  surface-to- 
volume  ratio  of  the  liquid,  thus  increasing  the  specific  rates  of 
mass,  momentum,  and  heat  transfer  and  the  vaporization 
rate.  Even  restricting  ourselves  to  injectors  used  in  diesel 
engines,  many  designs  exist  but  the  most  common  consists  of 
a  group  of  cylindrical  holes  l(X)-3(X)  pm  in  diameter.’  * 

When  a  liquid  is  forced  through  a  cylindrical  hole  into  a 
gas,  many  modes  of  breakup  are  observed.  In  the  one  relevant 
in  internal  combustion  engines,  no  outer  intact  length  is  seen 
and  the  jet  starts  diverging  at  the  nozzle  exit.  All  other 
parameters  being  the  same,  this  regime,  which  has  been  called 
the  atomization  regime,  is  reached  at  high  injection 
velocities— of  the  order  of  100  m/s  for  the  fuels  and  con¬ 
ditions  of  engine  applications.  ’ 

Because  of  their  practical  importance,  many  aspects  of 
atomizing  jets  have  been  studied  extensively.  In  the  I930’s 
significant  data  were  collected  on  such  global  quantities  as 
downstream  drop  sizes,  tip  penetration  rates,  and  average 
spray  angles.*'*  More  recent  efforts  have  attempted  to 
determine  the  structure  of  atomizing  jets.  For  example,  the 
outer  part  in  the  immediate  vicinity  of  the  nozzle  exit’  *  and 
the  inner  part  in  the  same  region  ’  are  being  studied  in  detail. 
Recently,  laser  techniques  have  been  used  to  measure  drop 
velocity,  drop  size,  and  liquid  volume  concentrations  at 
ambient  pressure  in  gas  turbine-type,  dilute  sprays  from  air 
blast  atomizers  in  which  more  than  90*l'o  of  the  light  is 
transmitted.' 

In  this  paper  we  report  drop  velocity  measurements  within 
nonvaporizing,  steady,  diesel-type,  fast,  dense  sprays  from 
single-hole  cylindrical  nozzles  into  compressed  nitrogen  at 
room  temperature  in  regions  so  close  to  the  injector  that  light 
transmissivity  was  as  low  as  2% . 
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Conditions,  Apparatus,  and  Procedure 

I  n-hexane  was  injected  into  quiescent  nitrogen  at  room 

temperature  but  at  such  pressures  that  the  ratio  of  the  gas 
density  to  the  liquid  density  was  0.02S6  and  0.0732.  Three 
I  single-hole  round  nozzles  with  sharp  inlet  and  outlet  and  two 

injection  pressures  1 1 .0  and  26.2  MPa)  were  used.  The 
I  effect  of  nozzle  geometry  was  explored  by  varying  the 

diameter,  127  and  76.2  fivn,  and  the  length-to-diameter  ratio. 
1  and  4.  These  parameters,  summarized  in  Table  I,  were 
chosen  because  earlier  work’  had  shown  the  gas-to-liquid 
density  ratio  and  the  nozzle  geometry  to  be  the  most  im- 
,  portant  variables  in  the  initial  formation  of  the  spray  and  the 

[  injection  velocity  is  important  in  the  subsequent  development 

and  propagation  of  jets.  Using  the  LDV  system  in  the  90-deg 
scatter  mode,  it  was  found  that  reliable  drop  velocity 
measurements  were  possible  at  locations  characterized  by  very 
high  drop  number  densities  (2l0"/m’)  and  high  velocity 
gradients  (up  to  Sm/s/mm).  Thus,  measurements  were  made 
across  the  jet  as  close  as  2.3  cm  from  the  nozzle  exit. 

The  experimental  apparatus  consisted  of  a  spray  chamber, 
a  liquid  pressurization  system,  a  nozzle  assembly,  and  LDV 
optics  and  instrumentation.  A  schematic  diagram  of  the  spray 
chamber  with  the  liquid  pressurization  system  and  the  nozzle 
assembly  is  shown  in  Fig.  I  and  details  are  available  in  Refs.  6 
and  10. 

The  spray  chamber  was  constructed  from  several  cylin¬ 
drical  steel  sections,  19-cm  i.d.  and  90  cm  in  total  length.  The 
window  section  has  four  quartz  windows  10  cm  in  diameter. 
The  liquid  pressurization  system  was  designed  to  maintain 
constant  pressures  up  to  207  MPa  during  injection.  The  liquid 
pressure  before  injection  was  measured  with  an  AMINC047- 
18340  gage  and  during  injection  with  a  Kistler  307A  trans¬ 
ducer  with  a  frequency  response  of  up  to  240  kHz,  a  Kistler 
S04  charge  amplifier,  and  a  Tektronix  7313  storage 
oscilloscope.  During  data  acquisition  the  pressure  changed  by 
less  than  l^o. 

The  single-hole  round  nozzles  were  drilled  directly  into  the 
nozzle  units,  whic  were  made  of  AISI  303  stainless  steel,  and 
examined  under  a  scanning  electron  microscope  to  assure  that 
the  desired  inlet  sharpness  was  obtained  and  not  altered 
during  the  tests.  Surface  roughness  was  less  than  5*^o  of  the 
diameter.  The  nozzle  unit  and  a  typical  spray  are  shown  in 
Fig.  2. 

The  liquid  injection  velocity  was  calculated  from  the 
measured  injection  duration,  the  area  of  the  nozzle,  and  the 
total  amount  of  liquid  injected  [ Pinj  =  volume/x(d/2)''Arl. 
This  volume  (mass)  mean  velocity  was  reproducible  to  within 
±l<^'o.  The  velocity  profile  within  the  nozzle  was  not 
measured  due  to  the  nozzle's  smallness  and  the  very  high 
velocity  of  the  liquid,  but  is  expected  to  have  been  flat,  except 
near  the  nozzle  walls,  because  the  upstream  liquid  was 
essentially  quiescent  and  the  lengih-to-diameier  ratio  of  the 
nozzle  was  much  too  small  for  transition  to  turbulence  away 
from  the  nozzle  walls.  Using  the  measured  l’,„,  and  assuming 
ideal  velocity  except  near  the  walls,  we  estimate  the  injection 
momentum  to  have  been  C,p,  P^,|»-(<//2)’  with  C,  = 
1.0±  10%. 

Two  LDV  systems  were  used.  Both  systems  employed  a 
Lexel  95-2  argon-ion  laser  operated  at  powers  between  0.1  and 
I  W  and  TSl  dual  beam  focusing  and  W-deg  scatter  collection 
optics.  The  two  LDV  systems  used  different  laser 
wavelengths,  514.5  and  488.0  nm:  fringe  spacing,  6.17  and 
2.45  pm;  and  focusing  optics  focal  lengths,  600  and  250  mm, 
respectively.  The  second  system  (Fig.  3)  was  superior  because 
of  its  capability  to  measure  flow-  reversal  with  the  TSl  9180 
frequency  shift  unit  and  its  smaller  probe  volume 
(0.2x0. 2x0. 2  mm’  vsO.2xO.2xO  7  mm').  When  using  the 
system  with  the  Bragg  cell  the  LDV  signal  was  electronically 
downmixed  to  0  Hz  for  the  axial  velocity  component 
measurements  in  the  core  region  where  the  drop  velocity 
fluctuations  were  small  compared  to  the  mean  velocity  and  to 
5  MHz  near  the  outer  edge  of  the  spray  and  for  all  the 


measurements  of  the  radial  velocity.  The  beam  crossing  angle 
was  measured  both  with  and  without  the  Bragg  cell  to  an 
accuracy  of  ±0.SV».  The  parallelness  of  the  two  beams  before 
the  focusing  lens  was  also  measured  and  found  to  be  within 
0.03  deg.  The  effect  of  this  on  the  beam  crossing  location  was 
estimated  to  be  negligible.  The  distance  from  the  beam  waist 
to  the  beam  crossing  was  also  calculated  and  found  to  be 
smaller  than  the  beam  waist  length.  Therefore,  it  was  assumed 
that  the  beams  crossed  at  their  beam  waist. 

To  allow  for  precise  and  repeatable  positioning  of  the  probe 
volume  within  the  spray,  the  entire  optical  system  was 
mounted  on  an  X.  Y,Z  tranversing  table  while  the  spray 
chamber  remained  fixed.  The  reported  axial  drop  velocities 
were  measured  along  (he  -Z  axis  and  the  radial  drop 
velocities  along  the  -  Y  axis  (see  coordinate  system  in  Fig.  3), 
At  each  axial  location,  first  the  axis  of  the  spray  was  located, 
based  on  the  symmetry  of  the  measured  velocity  profiles,  and 
the  radial  distances  were  referred  to  it.  A  minimum  of  five 
measurements  were  made  between  the  centerline  and  the  half- 


Fig.  I  Spray  apparatus  and  details  of  the  injeelion  cylinder:  I)  in¬ 
jection  cylinder,  2)  spray  chamber,  3)  window,  4)  gas  reservoir,  5) 
driver  gas  pressure  gage,  6)  liquid  reservoir,  7)  bleeding  valves,  81  lest 
liquid  pressure  gage,  9)  band  pump,  10)  lest  liquid  tank,  1 1 )  drain.  12) 
chamber  gas  pressure  gage,  13)  nitrogen  cylinder,  14)  regulators,  IS) 
solenoid  valves,  16)  valve  unit,  17)  upper  control  gas  chamber.  18) 
lower  control  gas  chamber,  19)  control  gas  ports,  20)  lest  liquid 
conduit,  21)  nozzle  piece,  22)  leak  off,  23)  nozzle,  24)  lest  liquid. 


Fig.  2  Nozzle  piece  and  l>  pica)  spray  (length  of  the  field  270if). 
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Table  1 

Spra)  conditions* 

Series 

p,. 

MPa 

MPa 

m/s 

Nozzle. 
tHtim)  l  /d 

X/d 

A 

1.48 

0.0256 

It.O 

127 

127-4 

600,800 

B 

4.24 

0.0732 

II.O 

127 

127-4 

300.400,500,600 

C 

4.24 

0.0732 

26.2 

194 

l27-» 

400,500.600 

D 

4.24 

0.0732 

n.o 

149 

76-» 

300.600.700.800 

E 

1.48 

0.0236 

It.O 

I2S 

76- 1 

300 

‘Liquid:  n-hexane.  p 
Room  lempcraiurc. 

*665  kg 

m’.  II,  • 

3.2k  I0*N 

$/m^.  -  1  *4x  10^ 

N/m.  Gas:  niirogen 

Table  2  Errors  in  LDV  measuremenis'* 

Source 

Type 

Estimated  error  and  correction 

LDV  optics 

Beam-crossing  angle 

0.5V, 

Fringe  spacing 

No  error  due  to  refraction  index 

Bragg  cell 

<15’* 

LDV 

Clock  counting 

<  I**!  based  on  an  intrinsic 

electronics 

accuracy  of  1  ns 

Noise 

<3Vc  with  SNR  =  10  and  Nf,  «4 

Spray 

Velocity  gradient 

s  2V,  for  the  mean  value 
s5^o  for  fluctuation  amplitude 
based  on  the  maximum  velocity 
gradient  of  5(m/'5)/mm 

Effect  of  drop  size 

No  sensible  effect  on  velocity 

distribution 

distribution  was  found  for 
X^iOOd 

Spray*io*spray 

variation 

Negligible 

Velocity  biasing 

Correction  applied  for 

RFA,  ,<30»'o 

Statistics 

Axial  component 

s  2^«  for  the  mean  value 

S  3^«  for  fluctuation  amplitude 
based  on  2000  data  and  50^0 
relative  fluctuation  amplitude 

Radial  component 

No  data  reported  for  the  mean 
value  due  to  large  error 

radius;  therefore,  the  uncertainty  in  the  positioning  at  an  axial 
location  of  300c/  is  no  more  than  ±8.3'l’(i  of  the  half-radius 
and  at  axial  positions  further  downstream  the  positioning 
error  is  an  even  smaller  percentage  of  the  half-radius. 

The  frequency  was  measured  with  a  TSl  1990  counter 
processor  interfaced  to  a  Hewlett-Packard  21  MX  minicom¬ 
puter.  The  number  of  fringes  set  on  the  counter  processor 
over  which  the  frequency  was  measured  was  eight  or  four. 
Such  a  small  number  of  fringes  was  used  primarily  to  increase 
the  data  rate  in  the  region  of  high  drop  number  density.  The 
use  of  fewer  fringes  also  reduces  the  trajectory  bias  that  can 
occur  when  a  large  number  of  fringes  are  used  which  favors 
drop  trajectories  normal  to  the  fringe  planes.’ 

The  same  initiation  signal  that  operated  the  valve  unit  (item 
16  in  Fig.  1)  was  used,  through  a  control  module,  to  enable 
the  counter  to  transmit  data.  The  control  module  was  also 
used  to  insert  a  prefixed  delay  between  spray  initiation  and  the 
start  of  data  transmission  and  to  program  the  data  acquisition 
period  for  each  test.  The  delay  was  necessary  to  allow  the 
spray  to  reach  a  steady-state  condition  as  determined  from  the 
time  history  of  the  liquid  pressure. 

All  the  measurements  were  made  in  steady  sprays.  Even 
though  the  data  were  taken  in  a  0.5- 1.0  s  period  to  prevent 
recirculation  inside  the  chamber  from  affecting  the  spray,  the 
velocity  of  these  sprays  was  so  high,  and  their  size  so  small, 
that  their  characteristic  times  are  much  shorter  than  1 .0  s.  The 
longest  time  would  be  the  convection  time  at  the  farthest  axial 
station  (800c/)  for  the  largest  of  the  nozzles  (127  /im)  and  the 
slowest  of  the  injection  velocities  (127  m^s).  This  time  is  less 
than  0.03  s.  As  previously  stated,  during  the  measurements 


liquid  and  gas  pressure  were  constant  to  better  than  1%,  Far 
downstream,  the  0.5-1. Os  window  was  long  enough  to  collect 
more  than  20(X)  data  points  during  one  injection,  but  going 
upsteam  the  rate  of  acquisition  of  acceptable  data  decreased 
and  up  to  30  injections  were  necessary  to  obtain  the  same  total 
number  of  velocity  data.  For  each  set  of  measurements,  the 
number  of  rejected  data,  based  on  the  criterion  that  the  data 
should  fall  within  3.5o  of  the  mean,  was  always  less  than  1  t'o 
of  the  total. 

Error  Analysis 

An  extensive  analysis  of  experimental  errors  was  carried 
out.  A  detailed  discussion  is  available  in  Ref.  10.  Sources  of 
error,  their  estimated  magnitudes,  and  possible  corrections 
are  summarized  in  Table  2.  Only  the  discussion  of  the  effects 
of  drop  size  distribution  and  velocity  biasing  is  given  here. 

No  direct  measurement  of  the  drop  size  distribution  was 
made  from  which  the  possible  correlation  between  drop  size 
and  velocity  could  be  determined.  However,  along  the  axis  it 
was  found  that  for  X/c/s:300  the  velocity  distribution  is  not 
affected  by  varying  the  laser  power  or  by  limiting  the  observed 
range  of  drop  sizes  to  which  the  counter  processor  responded 
(achieved  by  reducing  the  upper  limit  of  the  dc  component  of 
the  Doppler  signals  accepted  by  the  LDV  counter  processor). 
Figure  4  shows  the  effect  on  the  mean  velocity  and  the 
fluctuation  amplitude  at  X/d=8(X)  of  a  variation  of  the 
amplitude  limit  setting  on  the  counter  processor  from  100  to 
2.  No  significant  change  is  seen  in  either  L,  or  fluctuation 
amplitude  whereas  the  data  rate  changed  drastically,  in¬ 
dicating  a  large  selectivity  in  drop  size.  It  was  estimated  that 
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Fig.  3  LDV  setup. 


Fig.  4  \  trillion  of  mean  and  fluctuating  components  of  the  drop 
axial  velocii)  vs  amplitude  limit  setting  of  the  counter  processor. 


at  least  90'^'o  of  the  total  number  of  observable  drops  were 
rejected  with  the  lower  setting. However,  upstream  of 
A'/d=300,  the  mean  velocity  begins  to  show  sensitivity  to 
variations  of  the  laser  power  (equivalent  to  changing  the 
amplitude  limit  setting)  and  at  Xld- 200  change  in  y,  is  quite 
large,  as  shown  in  Fig.  5.  For  this  reason  only  the 
measurements  for  X/c/z300  are  reported  and  discussed. 
From  the  point  of  view  of  the  structure  of  these  sprays,  the 
fact  that,  for  X/dz300,  the  drop  velocity  is  independent  of 
the  drop  size  along  the  axis  implies  that  drop  and  gas 
velocities  there  are  the  same;  that  is,  the  equilibrium  limit  has 
been  reached.  For  X/d<300,  along  the  spray  axis  it  is  con¬ 
cluded  that  the  drop  and  gas  velocities  are  different. 

The  major  source  of  uncertainty  in  the  estimates  of  errors 
was  the  velocity  bias  effects.  A  velocity  bias  error  can  occur  in 
fluctuating  flows  since,  with  uniform  particle  concentration, 
more  particles  are  sampled  per  unit  time  when  the  gas  velocity 
is  higher.  The  magnitude  of  the  velocity  bias  error 
depends  on  the  local  fluctuation  amplitude.  Buchhave” 
calculated  the  velocity  bias  errors  of  the  mean  and  the  fluc¬ 
tuating  velocities  for  three-dimensional  Gaussian  isotropic 
turbulence,  for  relative  turbulence  intensities  up  to  80^o, 
including  the  effects  of  trajectory  bias  but  without  frequency 
shifting  and  showed  that  the  errors  increase  monotonically 
with  the  turbulence  intensity.  At  30'?'o  turbulence  intensity  the 
velocity  bias  errors  were  between  5  and  ICl^o  for  the  mean 
velocity,  and  between  3  and  7%  for  the  fluctuation  amplituuc, 
depending  on  the  ratio  of  the  minimum  number  of  fringe 
crossings  required  for  a  measurement  to  the  maximum 
number  of  fringe  crossings  available. 
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Fig.  5  Variation  of  mean  drop  axial  velocity  vs  laser  power. 
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The  determination  of  unbiased  averaged  quantities  from 
velocity  biased  data  can  be  achieved  only  if  the  time  each 
particle  remains  in  the  probe  volume  is  known  exactly.’  A 
simplified  one-dimensional  a  posteriori  correction  can  be 
applied  by  weighting  each  individual  measurement  with  a 
weighting  factor  inversely  proportional  to  the  measured 
velocity  component."  This  correction  method  can  reduce  the 
velocity  bias  error  when  the  relative  fluctuation  amplitude  is 
small.  For  example,  after  correction,  at  30%  turbulence 
intensity,  the  velocity  bias  errors  are  less  than  4.3%  for  the 
mean  velocity  and  between  2  and  S.S%  for  the  fluctuation 
amplitude."  For  our  sprays  the  one-dimensional  correction 
method  was  considered  acceptable  and  was  applied  as  long  as 
the  relative  fluctuation  of  the  drop  velocity  was  less  than  or 
equal  to  30%.  Above  30%,  the  residence  time  can  no  longer 
be  approximated  by  the  measured  velocity  component  and  the 
one-dimensional  correction  overcorrects  the  mean  velocity 
and  introduces  errors  that  can  be  larger  in  magnitude,  and 
opposite  in  sign,  than  those  of  the  uncorrected  data.  Thus,  the 
correction  does  not  improve  the  accuracy  of  the  data  and  the 
uncorrected  data  may  be  in  error  by  as  much  as  22%  for  the 
mean  velocity  and  10%  for  the  fluctuation  amplitude  at  80% 
turbulence  intensity  if  the  angular  dependence  is  minimized  by 
frequency  shifting." 

Besides  the  drop  velocity  distribution,  the  half-radius  as 
determined  from  the  drop  axial  mean  velocity  is  also  affected 
by  the  velocity  bias  error  because  at  the  half-radius  location 
the  fluctuation  amplitude  is  about  43%.  The  half-radius  was 
derived  both  from  data  which  had  been  velocity  bias  corrected 
around  the  half-radius  location  (fully  corrected)  and  from 
data  to  which  no  velocity  bias  correction  had  been  applied 
around  the  half-radius  location  (partially  corrected).  In  both 
cases  the  one-dimensional  velocity  bias  correction  was  applied 
for  relative  fluctuations  smaller  than  30%,  i.e.,  near  the  spray 
centerline.  The  two  sets  of  half-radius  results  differ  from  each 
other  by  10%.  Buchhave"  showed  that  without  correction  the 
half-radius  is  overestimated  and  with  correction  it  is  un¬ 
derestimated.  Thus,  based  on  half  the  difference  of  these 
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results  the  uncertainty  in  the  half-radius  due  to  the  velocity 
bias  error  is  about  ±  S^o.  But  when  all  other  sources  of  error 
are  considered,  the  total  uncertainty  in  the  half-radius 
becomes  ±  6.6*70 . 

The  velocity  bias  correction  is  strictly  valid  only  if  the  drop 
number  density  can  be  assumed  to  be  uncorrelated  with  the 
local  velocity  and  provided  no  other  limitations  on  the  mean 
sampling  rate  interfere  and  further  complicate  the  problem. 
The  drop  concentration  also  must  be  sufficiently  low  so  that 
most  of  the  time  there  is  only  one  drop  in  the  probe  volume. " 
This  condition  was  only  marginally  satisfied. 

In  summary,  the  one-dimensional  correction  was  applied 
on  the  drop  velocity  data  as  long  as  the  relative  fluctuation 
amplitude  was  less  than  or  equal  to  The  result  was  that 
the  total  maximum  error  of  the  reported  mean  velocity  was 
estimated  to  be  less  than  lO^o  up  to  the  half-radius  and  in¬ 
creasingly  larger  outward  where  it  can  be  as  large  as  40*7o  for  a 
relative  drop  velocity  fluctuation  amplitude  of  80Ve.  The  total 
error  of  the  reported  fluctuation  amplitudes  is,  in  general, 
smaller  than  that  of  the  mean  velocity. 

Results  and  Discussion 

The  parameters  varied  in  the  experiment  are  listed  in  Table 
I .  Five  conditions  were  examined  with  different  combinations 
of  two  gas-to-liquid  density  ratios,  two  injection  pressures, 
and  three  nozzle  geometries.  For  each  condition,  the  axial  and 
radial  components  of  the  drop  velocity  were  measured  by 
LDV  at  several  axial  and  radial  locations.  The  measurements 
discussed  in  this  section  were  taken  with  300-800  nozzle 
diameters  from  the  nozzle  exit  (2.29-10.16  cmj. 

We  will  consider,  first,  average  drop  velocities;  then, 
fluctuations  of  the  drop  velocities;  and,  finally,  the  im¬ 
plications  of  the  measurements  about  the  structure  of  these 
sprays.  In  the  figures,  the  corresponding  fluid  quantities 
measured  by  \\  ygnanski  and  Fiedler'^  by  hot-wire 
anemomeiry  in  isothermal,  low  Mach  number  (in¬ 
compressible)  air-into-air  jets  are  shown  for  reference  and 
comparison,  but  without  implying  that  they  constitute 
definitive  measurements  of  incompressible  jets. 

Figure  6  shows  that,  at  sufficient  distance  from  the  injector, 
the  ratio  of  the  centerline  velocity  to  the  injection  velocity 
tends  to  depend  only  on  Xld  .  In  this  figure,  data 

other  than  ours  are  from  Ref.  15.  According  to  Hinze,'*  this 
scaling  was  first  proposed  by  Thring  and  Newby,'''  but  was 
also  obtained  analytically  by  Kleinstein.^° 


X  /  w(a,/^)-*i 

Fig.  6  Cenlerlinf  vrlneity  dtci)  of  different  jets  ts  dimensionless 
axial  distance. 


To  understand  the  possible  reasons  for  the  scatter  of  the 
data  in  Fig.  6,  consider  a  steady  injection  of  a  fluid  f  into  a 
second  fluid  g.  At  sufficient  distance  from  the  injector  the 
flow  rate  of  the  injected  mass  becomes  negligible  with  respect 
to  that  of  the  entrained  mass  and  the  subsequent  development 
of  the  jet  depends  on  the  ambient  fluid  entraining  more 
ambient  fluid  as  in  “incompressible”  jets.  This  fully 
developed  incompressible  jet  limit  must  eventually  be  reached 
independently  of  the  nature  of  the  fluids  and  of  the  structure 
of  the  development  region.  Then,  neglecting  viscous  stresses 
and  pressure  gradients,  conservation  of  axial  momentum 
gives 


C,p,P'S,cirf'/‘'  =  C,P,U{.„rj5., 


C  ydjpjpX' 

Vo,p~^4C,C^ )  X-Xo 

where  C,  relates  the  self-preserving  distributions  of  the  mean 
axial  velocity  and  velocity  fluctuations  to  the  mean  centerline 
velocity  and  is  equal’^''*  ”  to  0.846  ±2.9^7o;C= 0.0868  ±7.9% 
relates  the  half-radius  to  the  distance  from  Xq,  the  virtual 
origin'*-^''^*;  and  C,  relates  the  injection  momentum  to  the 
injection  centerline  velocity.  For  uniform  injection  velocity 
profile  C,  =  I  and  the  coefficient  of  Eq.  (2)  is  6.3  ±9%  (Capp 


(!) 

(2) 


Fig.  7  Dimensionless  hiir-nidius  of  the  drop  axial  velocity  vs 
dimensionless  axial  distance. 


Fig  a  Sell-xtmllar  prnnie  nf  the  dimenxinniexx  mean  drop  axial 
vrlocti^ 
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and  George.^’  for  example,  give  5.8).  In  Fig.  6  the  continuous 
line  corresponds  to  A'(,  =  0  whereas  the  two  dotted  lines  cor¬ 
respond  to  {X„l<l)/{p,/p,)  =  30/(13. 7)  -  and  1(X)/{39.1)'' 

=  8.12  and  16.0,  respective!) .  The  Xo  values  of  30  and  100 
will  be  explained  later.  The  three  lines  show  that,  at  sufficient 
distance  from  the  injector,  the  structure  and  length  of  the 
development  region  become  immaterial. 

Often  there  is  uncertainty  about  the  value  of  C,  for  dif¬ 
ferent  experiments  because  the  injection  momentum  generally 
is  not  measured  and  because  different  authors  may  use 
velocities  other  than  the  centerline  velocity  as  reference  in¬ 
jection  velocity.  As  previously  stated,  we  did  not  measure  the 
injection  velocity  profile.  We  measured  the  mass  mean  in¬ 
jection  velocity,  and  estimated  that  the  injection 

momentum  was  equal  to  C,p,  F‘’„j-ird^/4  with  C,  =  l±10"?o. 
Thus,  using  for  in  Fig.  6  far  from  the  injector  our 
data  must  fall  on  the  continuous  line.  For  data  other  than 
ours  we  used  the  injection  velocity  given  by  the  various 
authors.  We  believe  that  the  main  reason  for  the  scatter  of  the 
data  at  large  values  of  X/d  {p,/p, ) is  the  uncertainty  about 
the  injection  momentum,  i.e.,  about  the  value  of  C,. 

Another  factor  could  be  that  the  fully  developed  local 
equilibrium  condition  was  not  reached  by  all  the  jets  of  the 
figure,  as  suggested  by  Shearer  and  Faeth.”  For  A'/rfa300, 
earlier  it  was  shown  that  drop  and  gas  velocity  tend  to  be  the 
same  along  the  axis  and,  presently,  it  will  be  shown  that  mean 
and  fluctuating  components  of  the  drop  velocity  have  reached 
self-preserving  profiles.  Thus,  the  conditions  for  convergence 
to  the  solid  line  of  Fig.  6  are  met  and  our  data  tend  to  fall  on 
it. 


Finally,  accurate  measurements  in  these  jets  continue  to  be 
difficult  and  experimental  errors  may  be  present.  For 
example,  the  data  of  Wygnanski  and  Fiedler  are  seen  to  tend 
to  a  different  slope  in  Fig.  6.  but  Capp  and  George”  report 
that  their  faster  decay  may  be  due  to  wall  effects. 

Equation  (2)  can  now  be  used  to  determine  somewhat  more 
precisely  the  virtual  origin.  If  Fm/U.  ti  is  plotted  vs  X  in 
linear  scales,  our  far-field  data  must  converge  on  straight  lines 
of  known  slope,  corresponding  to  the  continuous  line  of  Fig. 
6.  By  extrapolating  backward,  the  intersections  of  these  lines 
determine  the  virtual  origins.  We  found  Xo  =  iOd  for 
p,/p,  =  13.7  and  Xg<=\00d  for  p,/pj  =  39.1,  the  other 
variables  having  no  clear  effect  on  Xg  within  the  accuracy  of 
our  data.  Obviously,  the  convergence  to  fully  developed  jets 
occurs  further  downstream  in  our  sprays  than  in  in¬ 
compressible  jets  {p,/p,  =  1.0)  for  which  the  reported  Xg/d  is 
always  smaller  than  10.  Having  the  virtual  origins,  the  value 
of  C  was  determined  by  curve  fitting  the  measured  half-radius 
data  (Fig.  7),  and  found  to  agree  with  the  value  obtained  from 
incompressible  jet  data. 

Another  implication  of  Eq.  (2)  is  that  the  tip  penetration 
rate  of  sprays  for  constant  pressure  injection  into  a  stationary 
unconfined  gas  can  be  estimated  using  the  injection 
momentum  and  the  nozzle  diameter,  and  without  detailed 
knowledge  of  their  structure.  The  estimate  is  only  rough  near 
the  injector  but  becomes  accurate  far  from  it.  This  is  because 
under  such  conditions  the  velocity  of  the  tip  of  a  transient 
spray  is  close  to  70%  of  the  local  steady-state  centerline 
velocity.”'” 
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So  far  we  have  considered  only  the  mean  centerline  velocity 
and  the  half-radius.  Figure  8  shows  that  for  A'/c/>300  the 
parallel  between  the  mean  axial  drop  velocity  in  our  sprays 
and  the  mean  axial  fluid  velocity  of  incompressible  jets  ex¬ 
tends  to  all  axial  and  radial  locations  and  all  conditions. 


c)  Platnrss. 

Fig.  10  Seir-similir  profiles  of  the  dimensionless  fluctuillon  of  the 
drop  radial  telocils. 


The  parallel  extends  also  to  the  fluctuating  components  of 
drop  and  fluid  velocities  but  with  increasing  uncertainty.  At 
each  axial  and  radial  location  the  measured  axial  and  radial 
components  of  the  drop  velocity  show  a  distribution  of 
values.  The  standard  deviation  and  skewness  and  flatness  of 
these  distributions  are  shown  in  Figs.  9  and  10.  It  is  seen  that 
near  the  axis,  the  amplitude  of  the  velocity  fluctuation  of  the 
drops  falls  into  the  lower  part  of  the  range  of  incompressible 
jets  and  appears  to  be  larger  at  the  edge.  Near  the  axis  the 
distribution  of  the  drop  velocity  is  also  Gaussian,  as  indicated 
by  the  skewness  and  Harness.  At  the  edge  of  the  sprays  the 
shapes  of  the  spray  drop  and  incompressible  jet  velocity 
distributions  appear  to  differ.  Those  of  the  drops  tend  to 
remain  more  Gaussian. 

The  effects  of  corrections  and  differences  between  ex¬ 
perimental  techniques  on  the  disparity  between  the  drop 
velocities  at  the  edge  of  our  sprays  and  the  Huid  velocity  at  the 
edge  of  incompressible  jets  were  evaluated  for  the  case 
A^  =  26.2  MPa, />,  =4.24  MPa,  nozzle  127-4,  A'/rf=400,  and 
r=-l.lSros-  Results  are  shown  in  Fig.  9  (numbered  1-4  at 
r/To  5  =  2.75)  and  Table  3.  The  measured  drop  velocity 
distribution  is  nearly  Gaussian  without  the  one-dimensional 
correction  (results  I  in  the  skewness  and  flatness  graphs)  but  if 
the  one-dimensional  correction  is  applied  it  becomes  much 
narrower  (test  2  in  Table  3)  and  its  parameters  move  closer  to 
the  incompressible  jet  values  (results  2  in  Hatness  graph). 
Since  the  real  quantities  are  most  likely  between  the  corrected 
and  the  uncorrected  ones,  there  is  a  tendency  toward  better 
agreement. 

The  inHuence  of  the  measuring  technique  appears  to  be 
even  stronger  and  to  tend  to  further  close  the  gap  between  our 
results  and  those  of  the  incompressible  jet.  It  is  recalled  that 
most  of  the  incompressible  jet  data,  including  those  of 
Wygnanski  and  Fiedler,  were  taken  with  hot-wire 
anemometry.  We  can  approximate  the  spray  data  that  would 
have  been  obtained  by  hot-wire  anemometry  by  disregarding 
the  direction  of  the  drop  velocity  as  determined  by  the  Bragg 
cell.  The  spray  data  thus  obtained  would  move  closer  to  the 
corresponding  incompressible  jet  data  as  shown  by  results  3 
and  4  in  Table  3  and  Fig.  9.  Intermittency  can  also  contribute 
to  explaining  the  differences.  Our  drops  come  to  the  edge  of 
the  spray  primarily  with  the  gas  from  the  core  and  reHect  the 
Huctuations  of  this  gas.  Fluctuations  that  are  averaged  also 
over  the  potential  Huid  would  be  of  smaller  amplitude. 

Thus,  it  is  likely  that  the  difference  between  our  drop 
velocity  at  the  edge  of  our  sprays  and  the  corresponding  fluid 
velocity  at  the  edge  of  incompressible  jets  is  less  than  show'n  in 
Fig.  9.  However,  correct  values  cannot  be  established  at  this 
time. 

In  summary,  indications  are  that  for  X/d>300,  drops  and 
gas  are  nearly  in  equilibrium.  As  a  further  elaboration  of  the 
question  of  equilibrium,  estimates  of  the  characteristic 
relaxation  times  of  drop  velocity^  and  gas  turbulence*’  were 
considered.'®  It  was  concluded  that  at  A’=300</the  condition 
fot  local  equilibrium  should  be  met  at  least  by  the  more 
numerous,  smaller  drops.  Downstream  of  A'=  300tf,  the  drop 
relaxation  time  does  not  increase  appreciably,  because 
collisions  and  coalescence  are  no  longer  numerous  at  low 
liquid-volume  fractions,^*  whereas  the  turbulent  eddy  time 
continues  to  increase  as  XK  Thus,  a  station  is  reached  at 


Table  3  Effects  of  the  Bragg  cell  and  the  one-dimensional  correction  on  the  drop 
velocll)  dislributinn 


Tesi 

Brapg  cell 

C  orrcciion 

m/s 

m  s 

s. . 

F,., 

1 

Yes 

No 

0  024 

o.sot 

0.57 

5  32 

2 

Yes 

ID 

0  003 

0  I9R 

0.63 

18  44 

3 

Nl3 

No 

0  360 

0.351 

1.89 

7  63 

4 

No 

ID 

0  109 

0  166 

3.82 

25.16 
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which  most  of  the  drops  are  in  equilibrium  with  most  of  the 
eddies. 

Finally,  it  is  recalled  that  for  X<300d  all  indications  are 
that  local  equilibrium  is  not  present.  This  lack  of  equilibrium 
near  the  injector  and  a  successive  but  selective  equilibration 
process  downstream  are  in  general  agreement  with  those 
described  by  Faeth.” 

Summary  and  Conclusions 

Axial  and  radial  components  of  the  drop  velocity  were 
measured  by  LDV  at  numerous  radial  and  axial  locations 
(300sA'/ds800)  within  diesel-type  n-hexane  sprays  but  at 
room  temperature,  under  steady  injection  pressure,  and  into 
initially  quiescent  nitrogen.  Five  configurations  were  em¬ 
ployed  that  differed  in  gas-io-liquid  density  ratio  (0.0256, 
0.0732),  injection  velocity  (127, 149,  and  194  m/s),  and  nozzle 
geometry  (straight  cylindrical  holes  with  sharp  inlet  and 
outlet,  diameters  of  76  and  127  ^m,  and  length-to-diameter 
ratio  of  4  and  1 ).  The  LDV  system  was  made  up  of  an  argon- 
ion  laser,  TSl  dual  beam  LDV  optics  with  a  Bragg  cell,  90-deg 
scatter  collection  optics,  counter  processor,  and  minicom¬ 
puter. 

In  the  error  analysis  we  considered  errors  related  to  the 
LDV  optics,  LDV  electronics,  velocity  gradient  broadening, 
drop  size  distribution,  drop  number  density,  probe  volume 
position,  finite  sample  size,  and  velocity  bias.  The  magnitude 
of  the  velocity  bias  error  was  particularly  difficult  to  estimate 
but  it  was  concluded  that  for  the  reported  data  the  total  error 
is  less  than  lO^o  from  the  centerline  to  the  half-radius  and 
progressively  greater  and  more  difficult  to  evaluate  beyond 
the  half-radius.  Significant  uncertainties  about  the  magnitude 
of  errors  at  the  outer  edge  of  jets  w  ith  either  LDV  or  hot-wire 
anemometry  suggest  that  this  region  is  still  poorly  charac¬ 
terized  even  for  incompressible  jets. 

It  was  found  that  at  sufficient  distance  from  the  nozzle  so 
much  gas  has  been  entrained  that  the  subsequent  structure  of 
the  spray  is  dominated  by  the  entrained  gas  and  that  the  fully 
developed  incompressible  jet  structure  and  the  equilibrium 
limit  are  reached.  This  condition  is  being  approached  3(X) 
nozzle  diameters  from  the  nozzle  as  shown  by  all  the 
measured  drop  velocity  parameters:  jet  half-radius,  centerline 
velocity  decay,  axial  mean  velocity  distribution,  axial  and 
radial  velocity  fluctuation  distributions,  and  independence  of 
drop  velocity  on  drop  size.  However,  the  condition  is  not 
verified  at  distances  smaller  than  300  nozzle  diameters. 

For  engine  applications,  the  distance  up  to  300  nozzle 
diameters,  i.e.,  the  first  3-12  cm,  is  the  one  of  prime  interest. 
If  it  is  true  that  vaporization  should  speed  up  the  achievement 
of  the  fully  developed  state,  it  is  also  true  that  neither  the 
ambient  gas  nor  the  injection  pressure  are  steady  during  fuel 
injection  and,  therefore,  accurate  computations  of  engine 
sprays  are  likely  to  require  knowledge  of  the  non-equilibrium 
developing  region. 

On  the  other  hand,  under  steady  conditions  it  is  possible  to 
estimate  the  steady-state  centerline  velocity  and  the  tip 
penetration  rate  for  a  broad  range  of  density  ratios  using  only 
the  injection  momentum  and  the  nozzle  diameter.  The 
estimate  becomes  accurate  far  from  the  injector. 
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ABSTRACT 


A  better  understanding  and  characterization  of  the  formation  and  propaga¬ 
tion  of  high  velocity  sprays  from  single  hole  cylindrical  nozzles  is  of  impor¬ 
tance  both  fundamentally  and  practically.  The  steady  and  transient  structure  of 
these  sprays  is  qualitatively  similar  to  that  of  Incompressible  Jets  but  the 
breakup  of  the  liquid  column  into  drops  and  the  presence  of  drops  introduce 
substantial  quantitative  differences*  Measurements  of  the  angle  of  the  spray 
and  of  the  size  of  the  drops  near  the  nozzle  suggest  that  the  breakup  of  the 
outer  surface  of  the  liquid  Jet  is  due  to  aerodynamic  forces  that  lead  to  the 
rapid  and  selective  growth  of  surface  perturbations  generated  within  the  nozzle. 
The  state  and  mechanism  of  disruption  of  the  Inner  part  of  the  liquid  Jet  Is 
less  clear  but  sufficiently  downstream  only  individual  drops  are  present. 

Recent  LDV  drop  velocity  measurements  and  detailed  multidimensional  computations 
have  shown  that  at  distances  of  the  order  of  hundreds  of  nozzle  diameters  so 
much  ambient  gas  has  been  entrained  by  the  spray  that  the  subsequent  structure 
of  the  Jet  is  dominated  by  the  entrained  ambient  gas  and  the  fully  developed 
incompressible  Jet  structure  and  drop-gas  equilibrium  are  approached. 
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NOMENCLATURE 


A,  B,  C,  Cjj,  Cn,  Ct 


Dimensionless  constants  of  Eqs.  15,  16,  17, 
1,  2,  3,  respectively 
Nozzle  discharge  coefficient 


Nozzle  diameter 


Drop  diameter 
Turbulent  dlffuslvlty 


Flatness  of  the  fluctuation  of  the  axial 


component  of  the  drop  velocity,  v^*'*/ Vj{'‘ 
Nozzle  length 

Specific  momentum  flow  rate 


Pressure 


Volumetric  flow  rate 

Jet  half  radius  (half  the  width  at  half  the  depth) 

Liquid  jet  Reynolds  number,  d/vi 

Skewness  of  the  fluctuation  of  the  axial 

_ ^ _ _  3/2 

component  of  the  drop  velocity, 

Gas  velocity  vector  of  components  u^i 

Mass  mean  Injection  velocity 

Drop  velocity  vector  of  component  Ux, 

Jet  Weber  number,  p£  d/o 


Axial  distance  from  the  nozzle  exit 

Virtual  origin 

End  of  Intact  liquid  core 

Beginning  of  fully  developed  entrained  gas  jet 
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SUPERSCRIPTS 


* 

SUBSCRIPTS 

CL 
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Radial  distance 

Injection  pressure  difference 

Spray  angle 

Angle  corresponding  to  tQ.s 
Gas  viscosity 
Liquid  viscosity 
Surface  tension 


Standard  deviation  of  the  fluctuation  of  the  axial 

1/2 


component  of  the  gas  velocity,  Ux'^ 


Standard  deviation  of  the  fluctuaton  of  the  axial 

_ -  1/2 

component  of  the  drop  velocity,  Vx''^ 


Gas  density 
Liquid  density 
Characteristic  time 


Fluctuation 

Applicable  to  far  field  of  sprays 

On  the  centerline,  r  ■  0 
Entrained 


0 


At  the  nozzle  exit,  x  -  0 
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INTRODUCTION 

Sprays  exhibit  a  large  variety  of  geometrical,  dynamic,  and  thermodynamic 
configurations  The  general  features  of  these  configurations  are 

understood  but  the  details  are  not  and  predictions  remain  elusive.  In  this 
field,  as  in  many  others,  the  stability  of  complex  Interfaces,  non-equlllbrium 
thermodynamics.  Intricate  reactions,  and  turbulence  ultimately  set  the  limits  of 
our  knowledge. 

This  Is  true  even  for  the  simplest  of  configurations,  that  Is  the  one  of 
Interest  In  this  review.  In  which  a  liquid  Is  Injected  Into  a  gas  through  a 

n 

single  straight  hole  of  circular  crosssection,  the  two  media  have  negleglble 
angular  momenta,  and  their  thermodynamic  states  are  such  that  vaporization  and 
chemical  reactions  can  be  Ignored.  Only  dynamic  forces  control  the  field  and 
the  evolution  of  Interfaces  and  turbulence  present  the  greatest  difficulties. 

A  general  Idea  of  the  structure  of  our  sprays  can  be  obtained  by  con¬ 
sidering  the  similar  and  better  understood  family  of  Incompressible  Jets 

Figure  1  shows  the  Initial  propagation  of  a  turbulent  Incompressible  Jet 
and  of  a  spray  The  propagation  of  the  Incompressible  Jet  Is  marked  by 

''^the  advancement  and  growth  of  a  head  vortex  that  Is  fed,  from  Its  downstream 
side,  by  the  Injected  medium  and  the  ambient  medium  that  was  entrained  in  the 
region  between  the  vortex  head  and  the  nozzle.  Practically,  this  Intermediate 
region  Is  In  Its  steady  state  configuration.  That  is,  the  head  vortex  leaves  a 
turbulent,  steady.  Incompressible  Jet  behind  Itself.  In  the  laboratory  frame, 
the  vortex  moves  at  a  fraction  of  the  local  steady  state  centerline  velocity. 
Within  the  steady  part,  the  shear  layer,  the  potential  core,  and  Its  end  can  be 


seen. 
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The  corresponding  picture  of  a  spray  raises  several  questions.  Its  Inter¬ 
nal  structure  cannot  be  seen  and  Is  not  known  precisely  because  no  technique  has 
been  found  to  probe  It  without  altering  It  significantly.  There  Is  a  head 
structure  but  Its  details  are  not  obvious.  It  Is  expected  that  when  the  head 
structure  Is  at  sufficient  distance  from  the  nozzle,  enough  ambient  gas  has  been 
entrained  to  set  up  a  flow  field  similar  to  that  of  incompressible  Jets  behind 
itself.  But  closer  to  the  nozzle,  the  liquid  core,  liquid  ligaments,  and  drops 
impose  different  configurations  on  the  shear  layer  and  on  the  head  of  the  Jet 
due  to  their  different  modes  of  momentum  transfer.  The  length  of  the  transition 
region  Is  not  known.  Behind  the  head,  however,  the  spray  divergence  angle, 
which  Is  the  only  quantity  that  can  be  measured  with  relative  ease,  rapidly 
achieves  its  steady  state  value  This  indicates  that  also  In  sprays  the 

adjustment  to  steady  state  occurs  primarily  within  their  head  region.  In  the 
steady  part  of  the  spray  of  Fig.  1,  as  In  the  corresponding  Incompressible  Jet, 
we  expect  the  mean  axial  velocity  to  point  downstream  and  the  mean  radial  velo¬ 
city  to  point  toward  the  axis  always.  Thus  we  call  them  full-cone  sprays  in 
contrast  to  hollow-cone  sprays  that  exhibit  mean  recirculation  flows  along  the 
ax! s  and  correspondingly  larger  spray  angles. 

Ir  Fig.  lb  the  spray  appears  to  start  diverging  Immediately  at  the  nozzle 
exit.  Although  this  Is  the  mode  of  breakup  of  interest  in  this  review,  there 
are  oth^r  modes  and  some  are  shown  In  Fig.  2.  All  other  parameters  being 
the  sam(,  at  very  low  Injection  velocities  a  Jet  breaks  up  many  diameters 
downstreams  and  forms  large  drops.  Surface  tension  is  the  disruptive  force 
(Fig.  2a).  As  the  velocity  Is  increased,  the  displacement  of  the  gas  by  the 
moving  undulated  liquid  surface  generates  a  pressure  distribution  tliat  aids  sur- 
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face  tension  in  amplifying  surface  waves  (Fig.  2b).  This  process,  which  is 
called  aerodynamic  Interaction,  eventually  becomes  the  main  dlstabillzing  force 
and  leads  to  the  formation  of  very  small  drops  while  being  opposed  by  surface 
tension  (Fig.  2c).  All  along  an  Intact  surface  is  visible  and  the  growth  of 
unstable  surface  waves  is  detectable.  These  are  classical  regimes  of  breakup  of 
circular  liquid  columns  and  major  contributions  to  our  understanding  of  them 
were  made  by  Lord  Rayleigh,  C.  Weber,  and  G.I.  Taylor.  They  continue  to  be 
the  subject  of  current  research  but  are  not  the  one  of  interest  to  us  here. 

The  regime  of  Interest  to  us  is  obtained  when  the  injection  velocity  is  so 
^^high  that  the  Intact  length  of  the  outer  Jet  surface  seems  to  disappear  and  the 
configuration  of  Figs.  2d  and  lb  is  obtained.  This  regime  is  called  the  atomi¬ 
zation  regime  and  the  forces  that  control  it  are  more  complex  and  less 

known  even  though  it  has  been  the  subject  of  extensive  research  because  of  its 
considerable  practical  Importance. 

Thus  we  will  concentrate  on  non-vaporizing,  non-reactive,  isothermal,  ato¬ 
mizing,  full-cone  sprays,  such  as  those  of  Figs,  lb  and  2d,  and  discuss  their 
possible  structure.  Ue  will  consider  first  the  steady  far  field  and  then  the 
breakup  and  development  regions.  For  liquids  such  as  water  and  hydrocarbon 
fuels  injected  at  room  temperature  into  standard  or  compressed  air  or  similar 
gases,  such  sprays  are  found  at  injection  velocities  of  the  order  of  10^  m/s. 

The  high  velocity  allows  large  mass  flow  rates  to  be  obtained  even  with  small 
nozzles.  Their  diameter  is  often  of  the  order  of  10^  pm  so  that  their  initial 
characteristic  time  is  of  the  order  of  10~^  s.  Drop  diameters  are  of  the  order 
of  10  gm.  The  combination  of  large  velocity  and  small  size  makes  detailed 
measurements  difficult  in  the  breakup  and  development  regions.  Thus  their  ini¬ 
tial  transient  was  observed  for  the  first  time  only  recently  (8). 
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THE  STEADY  FAR  FIELD 

The  classical  point  source,  boundary  layer,  similarity  solution  of  Tollmien 
for  incompressible  Jets  Is  a  satisfactory  guide  to  the  understanding  of  the 
global  structure  of  the  far  field  of  full  cone  sprays.  In  principle,  only  one 
constant,  such  as  that  relates  the  effective  turbulent  dlffuslvlty  D^  to  the 
specific  momentum  flow  rate  M,  is  left  indetermlned  by  this  solution  and  must  be 
evaluated  from  measurements.  Then  the  width  increases  as  x,  the  centerline 
velocity  decays  as  1/x,  and  the  ratio  of  the  entrained  mass  to  the  injected  mass 
increases  as  x: 


«  ■  V  “l.CL 

■io  ■  “O.CL 
“t  ■  =e  >V  “O.CL 

Tq  3  -  8  (x-Xq) 

=..CL  ■  »  C^)  d  :o.Cl/<’^»0> 

Q./Qo  ■  S 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


In  the  above  equations  the  constants  and  Cq  are  used  to  elate  the  spe¬ 
cific  momentum  flow  rate  and  the  volumetric  flow  rate  to  the  Injection  cen¬ 
terline  velocity  through  the  Injection  velocity  profile,  and  xq  Is  t’.e  virtual 
origin,  (In  practice  the  determination  of  Cgj,  xq  and  of  an  additional  needed 
constant  X2  is  nontrivial  even  for  Incompressible  Jets.  We  will  return  to  this 
subject  later.) 

But  at  this  point  the  growth  of  the  entrained  mess  and  local  characteristic 
convection  and  diffusion  times  are  of  Interest  to  uc.  Schilchting  gives 


*  O.Ui6i  and  Eq.  6  shows  that  for  C„  ■  Cq  1,  the  entrained  mass  is  already  10 


times  the  injected  one  at  22  nozzle  diameters  fror  the  virtual  origin.  Local 


convection  and  diffusion  times  can  be  obtained  dividing  (x-xq)  by  Ux, CL  fz'om  £q. 
5  and  r^o. 5  by  from  Eqs.  4  and  3*  The  two  times  are  of  the  same  order  and 
increase  as  (x-xq)^. 

For  our  sprays  the  injected  medium  is  a  liquid  and  the  entrained  medium  is 
a  gas,  so  that  the  initial  entrainment  process  can  be  expected  to  differ  from 
that  of  incompressible  Jets.  But  at  sufficient  distance  from  the  injector,  the 
Injected  mass  must  become  neglegible  In  comparison  to  the  entrained  one  because 
past  a  certain  distance  their  ratio  grows  as  x.  Downstream  of  such  region,  the 
momentum,  that  leads  to  the  entrainment  of  more  ambient  gas  resides  mostly  with 
/'^he  ambient  gas  that  was  entrained  earlier.  That  is,  the  structure  of  fully 
developed  incompressible  Jets  is  recovered. 

Imbedded  in  this  far  field  there  are  drops  that  move  in  equilibrium  with 
the  gas.  This  is  because  the  time  for  the  drop  velocity  to  relax  to  the  local 
gas  velocity,  Pid^^/lSpg^  has  an  upper  limit  whereas  the  fluid  times  continues 
to  grow  as  x^.  Since  there  is  a  distribution  of  drop  sizes  and  of  eddy  times, 
equilibration  will  be  selective  and  dependent  on  conditions  but  at  some 
appropriate  distance  from  the  injector,  all  sprays  become  incompressible  Jets 
^^ominated  by  the  entrained  ambient  medium  and  drops  move  within  them  as  markers 
of  the  motion  of  the  ambient  medium. 

Having  accepted  the  existence  of  this  limit,  Eqs.  1-6  can  be  modified  for 
direct  application  to  the  far  field  of  sprays.  Following  Thring  and  Newby 
and  Klelnstein  we  can  equate  the  axial  momentum  evaluated  at  the  nozzle 

exit,  where  the  density  is  P£,  to  that  evaluated  in  the  far  field,  where  the 
density  is  Pg  and  the  velocity  profile  is  fully  developed,  and  conclude  that,  as 
far  as  the  far  field  is  concerned,  a  liquid  spray  is  an  incompressible  Jet  with 
an  equivalent  specific  momentum  flow  rate  of  M(p£/pg).  Then  Eqs.  1-6  become 
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Q*/Qg  -  (Qjt/Qo>(pg/P4>^^^ 

d*  -  d(p^/Pg)^/2  (13) 

Parallel  Interpretations  are  that  in  the  far  field  a  spray  is  an  incompressible 
Jet  from  an  equivalent  nozzle  diameter  ^^2)  of  dCp^/pg)^/^  or  with  an  equivalent  ^ 
turbulent  dlffuslvlty  of  D4(p4/pg)l/2,  Notice  however  that  the  change  In 

the  virtual  origin  is  not  Identified  by  this  theory  since  this  quantity  is 
determined  by  the  structure  of  the  development  region.  The  ratio  of  the  tur¬ 
bulent  diffusion  time  to  the  drop  relaxation  time  becomes 


drop  i  t 

Several  authors  have  realized  the  existence  of  this  limit  and  taken  good 


advantage  of  it  But  only  recently  detailed  drop  velocity  measurements 

vlt'iln  sprays  of  this  family  have  provided  direct  evidence  of  it.  We  shall 
review  these  measurements  briefly,  but  first  a  note  of  caution.  In  specific 
applications,  particularly  to  time  varying  and/or  closed  volumes,  the  develop¬ 
ment  region  is  often  the  one  of  prime  interest.  Indeed  in  many  cases  no  more 
volume  or  time  are  provided  than  those  necessary  to  complete  some  specific 
degree  of  mixing  of  the  Injected  and  ambient  media.  In  the  far  field  limit,  the 
injected  medium  is  but  a  trace  within  the  ambient  one.  Thus  we  are  not  saying 


1 

or  Implying  that  the  local  equilibrium  limit  is  necessarily  achieved  in 
applicatlonsi  We  are  saying  that  these  limits  exist  and  are  conceptually  useful. 

The  data  of  Fig.  3  were  taken  in  the  sprays  of  Table  1  under  steady  injec¬ 
tions,  at  distances  greater  than  300d,  and  with  the  gas  contained  in  a  vessel  of 

/ 

such  size  that  a  wall-free  environment  was  approach  (1^).  These  sprays  looked 
very  much  like  those  of  Figs,  lb  and  2d.  Up  to  r  ■  tq, 5»  the  radial  profiles 
of  the  mean  axial  drop  velocity  and  of  the  amplitude,  skewness,  and  flatness  of 
the  fluctuation  of  the  drop  axial  velocity  are  seen  to  fall  within  the  range  of 
the  corresponding  fluid  quantities  measured  in  incompressible  jets  by  Wygnanski 
^^and  Fiedler  by  hot  wire  anemometry  and  by  others  with  hot 

wire  anemometry  and  laser  Doppler  veloclmetry.  The  same  parallel  was  also 
found  for  the  radial  component  of  the  drop  velocity,  for  the  centerline  velocity 
decay  and  for  the  width.  Also,  the  drop  velocity  was  measured  with  various 
laser  power  levels,  thus  weighing  the  measurements  in  favor  of  drops  of  dif¬ 
ferent  sizes,  and  the  results  were  found  to  be  Independent  of  it  for  x>300d.  If 
the  drops  had  not  been  in  equilibrium  with  the  gas,  their  velocity  would  have 
depended  on  their  size.  Thus  indications  are  that  for  the  full  cone  sprays  of 
^^^^able  1  the  ini.  ompressible  jet  and  drop-gas  equilibrium  limits  are  achieved 
around  x  ■  300d  and  for  r  <  ro.3.  For  spray  E,  at  x  ■  200d  the  measurements 
indicate  that  drops  are  not  in  equilibrium  with  the  gas  and  for  x>300d  but 
r>ro.5  the  evidence  is  inconclusive  because  of  large  errors  in  both  HWA  and  LDV 
data 

There  is  no  reason  to  expect  that  all  sprays  of  Table  1  should  reach 
equilibrium  at  the  same  axial  distance,  but  the  data  were  not  sufficiently 
numerous  or  sufflclenty  accurate  to  dlf ferenclate  and  the  300d  location  should 
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be  considered  Indicative.  By  extrapolating  backward  the  centerline  velocity 
decay,  and  also  checking  the  growth  of  the  Jet  width,  it  was  determined  that  xq 
■  30d  for  Pji/Pg  ■  13.7  and  xq  ■  lOOd  for  P£/pg  “  39.1.  If  the  300d  range  and 
the  two  values  of  xq  are  compared  with  the  corresponding  50d-70d  and  5d-10d 
values  for  Incompressible  Jets  the  reasonable  conclusion  is  reached  that 

full-cone  sprays  develop  into  incompressible  Jets  but  require  a  longer  distance 
and  that  such  distance  is  likely  to  be  an  increasing  function  of  p^/pg.  xhe 

density  ratio,  Pj^/pg,  is  seen  to  be  the  main  additional  parameter  for  the 
achievement  and  structure  of  the  far  field  of  these  sprays. 

THE  STEADY  NEAR  FIELD 

In  considering  the  near  field  we  will  start  from  the  outer  part  of  the 
spray  in  the  Immediate  vicinity  of  the  nozzle  exit.  It  is  of  Interest  to  know 
what  forces  breakup  the  liquid  surface  in  this  region.  From  photographs,  such 
as  Figs,  lb  and  2d,  it  appears  that  the  Jet  starts  diverging  immediately  at  the 
nozzle  exit.  Higher  resolution  images,  such  as  Fig.  6  which  will  be  recon¬ 
sidered  later,  show  Isolated  drops  and  an  opaque,  highly-lrregular,  diverging 
fluid  that  could  be  made  up  of  any  combination  of  drops,  ligaments,  blobs,  and 
deformed  intermingled  gas-liquid  continua. 

Hypothesis  about  the  breakup  mechanism  are  not  lacking  but  quan¬ 

titative  evidence  from  controlled  well  documented  unequivocal  experiments  has 
been  very  scant.  One  difficulty  is  that  too  many  events  that  can  contribute  are 
present  simultaneously  and  cannot  be  investigated  separately  while  still  con¬ 
sidering  the  same  family  of  sprays. 

For  example,  turbulence  of  the  liquid  was  suggested  as  the  main  dlstabl- 
llzlng  agent.  The  Reynolds  number  based  on  the  nozzle  diameter  is  generally 
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greater  than  10^  but  the  nozzle  length-to-diameter  ratio  is  seldom  greater  than 
6  and  fully  developed  turbulent  pipe  flows  are  not  present.  Moreover  increasing 
L/d  from  10  to  80  leads  to  smaller  divergence  angles,  l.e.  to  more  stable  Jets 
(19,20)^  Within  the  nozzle,  there  are  turbulent  wall  boundary  layers  and  the 
sudden  change  of  forces  within  them  that  occurs  at  the  nozzle  exit  has  been 
suggested  as  dlstablllzlng,  but  even  in  the  most  carefully  machined  metal 
nozzles  the  surface  roughness  Is  still  no  less  than  10  pm  with  diameters  In  the 
range  of  70  pm  to  300  pm.  It  Is  not  clear  whether  the  tall  ridges  and  deep 
valleys  formed  by  surface  roughness  trip  the  flow  or  trap  it.  In  any  case, 
^/^radluslng  the  exit  edge  with  various  curvatures  brought  about  no  measurable 

change  in  the  Initial  spray  angle  (20),  Cavitation  is  invariably  present  at  the 
entrance  of  practical  nozzle  of  this  family  and  has  been  suggested  as  the  main 
dlstablllzlng  agent.  But  cavitation  free  nozzles  have  been  found  to  give  Imme¬ 
diately  diverging  Jets  too  (20).  Their  angle  was  smaller,  indicating  greater 
stability,  but  the  geometry  of  the  cavitatlon-free  nozzles  was  so  different  from 
the  sharp-inlet,  sharp-outlet,  straight-wall  geometry  of  standard  nozzles  that 
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the  entire  nozzle  flow  field  was  also  different. 

Also  suggested  as  possibly  being  responsible  for  the  breakup  are  the 
rearrangement  of  the  cross-section  axial  velocity  profile,  and  liquid  supply 
pressure  oscillations.  But  none  of  these  mechanisms  alone  was  found  adequate  to 
explain  the  trends  exhibited  by  the  Initial  spray  angle  (^^»  20),  However,  any 
of  them  could  be  contributing  to  the  break  up  process  as  explained  below. 

Since  aerodynamic  Interaction  Is  known  to  cause  the  breakup  of  Jets  at 
lower  speeds,  the  suggestion  that  It  may  continue  to  do  so  at  higher  speeds  too 
Is  a  natural  one.  According  to  this  view  the  length  of  Jet  surface  over  which 


unstable  surface  waves  grow  becomes  shorter  and  shorter  as  the  speed  Is 
increased,  due  to  the  faster  growth  rate  of  the  unstable  waves,  until  It  becomes 
of  the  order  of  a  few  microns  and  is  no  longer  detectable.  Then  the  Jet  appears 
to  diverge  Immediately  at  the  nozzle  exit.  This  first  order  linear  perturbation 
theory  leaves  a  parameter  unspecified  that  can  be  interpreted  as  the  Initial 
amplitude  of  the  surface  perturbations  when  the  liquid  first  enters  the  gas.  If 
this  parameter  is  allowed  to  vary  with  nozzle  geometry,  as  it  would  be  the  case 
if  different  nozzle  flows  establish  different  initial  perturbation  levels 
through  some  combination  of  the  previously  mentioned  processes,  then  a  supple- 
rr>ented  aerodynarlic  interaction  mechanism  results  that  seems  to  comply  In  a 
iitting  way  with  a  rather  large  set  of  experimental  Information.  Castleman 
was  among  the  early  supporters  of  this  view.  Ranz  produced  a  theoretical 

framework  for  it,  extending  the  work  of  Taylor  (23),  ^.nd  Reitz  (24,11)  per¬ 
formed  a  comprehensive  evaluation  of  most  major  proposals  and  sharpened  the 
focus  on  it. 

The  supplemented  aerodynamic  theory  leads  to  the  prediction  (22,24,25)  Qf 
r'le  initial  spray  angle,  the  Initial  average  drop  size,  and  the  length  of  the 
nt.j.t  core.  The  equations  are  particularly  simple  in  the  limit 

’i\aj2/og  VJej^  >  1  which  is  the  one  of  practical  Interest  in  raany  cases  of  ato- 
.'Ization.  The  angle  is  determined  by  combining  the  radial  velocity  of  the 
fastest  growing  of  the  unstable  surface  waves  with  the  axial  injection  velocity: 
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(15) 


where-  A  is  a  constant  whose  value  depends  on  the  nuzzle  geometry  and  must  be 
determlnt'i  experimentally.  The  luirlal  average  drjD  diameter  is  assumed  to  be 
rroport  ion.il  to  the  lengtli  of  the  most  unstable  wj’  ■! 


where  B  Is  a  constant  of  order  one  and  Independent  of  the  nozzle  geometry  (but 
dependent  on  the  reference  velocity  that  is  chosen  for  the  break  up  process). 
And  the  length  of  the  intact  core  is  obtained  by  subtracting  the  mass  of  the 
drops  from  the  intact  liquid  column  as  they  are  formed: 


x.  p-  1/2 

“  C  (~) 


(17) 


where  C  is  a  proportionality  constant. 


The  prediction  that  in  the  atomization  regime  the  initial  value  of  the 
spray  angle  depends  almost  exclusively  on  the  density  ratio  and  the  nozzle 
geometry  is  surprising,  considering  the  many  parameters  that  could  effect  it, 
but  is  generally  born  out  by  measurement.  Figure  4  shows  that,  for  a  given 
nozzle  geometry,  large  changes  in  liquid  properties,  injection  velocity,  and  gas 
pressure  bring  about  only  very  small  systematic  trends  and  that  the  density 
ratio  dominates  (2^).  Figure  5  gives  an  example  of  the  effect  of  nozzle 
geometry,  all  other  parameters  being  the  same  In  spite  of  the  success 

/'^here  are  limitations:  The  mechanism  by  which  the  nozzle  geometry  influences 
the  angle  is  not  known;  although  the  angle  is  very  reproducible  in  any  given 
experiment,  its  value  depends  somewhat  on  definitions  and  measuring  techniques; 
when  a  broad  range  of  injection  velocities  is  explored  a  mild  trend  is  detected 
at  lower  density  ratios  that  does  not  conform  with  the  expected  one  (^0). 

Recently  (26,27)  ^^e  diameter  of  drops  was  measured  at  the  edge  of  ato¬ 


mizing  sprays  in  the  immediate  vicinity  of  the  nozzle  exit  using  photographs 
such  as  that  of  Fig.  6.  The  average  diameter  is  given  in  Fig.  7  for  the  con- 


dltlons  of  Table  2.  When  compared  with  Eq.  16,  the  measured  values  exhibit  the 
correct  trends  with  respect  to  injection  velocity,  liquid  properties  and  nozzle 
geometry  (It  should  have  no  effect  and  doesn't)  but  the  Incorrect  one  with 
respect  to  gas  density.  Moreover  the  measured  diameters  are  about  a  factor  of  3 
greater  than  expected.  However  It  would  appear  that  a  reasonable  explanation 
for  the  disagreements  exists  (2^).  The  drops  that  could  be  and  were  measured 
are  at  the  outer  edge  of  the  spray  and  not  at  their  formation  sites  to  which  Eq. 
16  applies.  Due  to  the  high  drop  number  density  In  the  region.  It  Is  unreason¬ 
able  to  expect  that  each  drop  clears  the  congested  area  and  comes  out  without 
colliding  with  other  drops.  When  collisions  and  coalescence  are  considered, 
with  a  model  that  will  be  mentioned  presently,  the  discrepancies  tend  to  disap¬ 
pear.  It  Is  true,  however,  that  one  set  of  measurements  and  computations  cannot 
possibly  be  considered  sufficient  to  close  this  complex  subject. 

Finally  the  data  of  Hiroyasu  et  al  support  the  prediction  of  Eq.  17 
with  respect  to  Injection  velocity  and  gas  density.  However  the  technique  they 
used,  based  on  measurement  of  the  electrical  resistance  between  the  nozzle  and  a 
screen  that  could  be  moved  axially,  detects  any  continuous  liquid  connection 
between  the  nozzle  and  the  station  of  the  screen.  A  ligament  or  a  mixture  of 
liquid  and  gas  would  give  the  same  signal  as  a  solid  liquid  core  free  of  gas. 
Thus  some  uncertainty  about  the  structure  of  the  core  persists. 

Individually  any  one  of  the  quoted  experiments  may  not  provide  adequate 
support  for  the  supplemented  aerodynamic  theory  of  etomization.  But  together 
they  fora  a  rather  consistent  picture.  By  comparison,  the  support  for  alter¬ 
native  theories  is  very  meager. 


THE  TRANSIENT 


Because  of  its  Importance  in  Diesel  engines,  many  measurements  have  been 
reported  of  the  velocity  of  the  tip  of  atomizing  jets  (28)  ^nd  It  Is  generally 
agreed  that  It  Is  about  70%  of  the  steady  state  centerline  velocity.  But 
details  are  anything  but  clear  as  Indicated,  In  the  first  instance,  by  the  lack 
of  consistency  In  the  reported  curve  fits  for  the  tip  velocity.  In  fact,  the 
transient  of  round  Jets  has  received  very  little  detailed  attention. 

There  seems  to  have  been  only  one  study  of  the  transient  of  laminar 
incompressible  Jets  (29)  concludes  with  wrong  time  and  length  scales  for 

^^lt.(^®^  Recent  measurements  (^^^  of  transient  Incompressible  Jets  were  limited 
in  value  by  complicated  nozzle  conditions.  Recent  numerical  studies  (2^))  of  the 
transient  of  laminar,  turbulent  and  spray  Jets  gave  reasonably  complete  Infor- 
mation  about  laminar  Jets,  because  there  are  no  uncertain  physical  parameters 
for  such  Jets,  but  only  indicative  Information  about  turbulent  and  spray  Jets, 
because  there  are  too  many  uncertain  parameters  for  them. 

The  diffusive  nature  of  Jets  makes  the  identification  of  their  most 
advanced  position  a  matter  of  definition  and  a  function  of  the  experimental 
technique.  However  It  would  appear  that  all  Jets  scale  up  In  time  at  least 
approximately.  That  a  scaling  exists  Is  shown  by  pictures,  such  as  those  of 
Fig.  1,  and  by  the  many  correlations  derived  from  them.  That  the  scaling  of 
sprays  is  complex,  or  Just  an  approximation,  is  Indicated  by  the  small  probabl- 
llty  of  the  existence  of  exact  similarity  solutions  of  their  constitutive 
equations. 

Indeed  If  we  reconsider  Fig.  lb  and  If  we  define  the  arrival  of  the  spray 
as  the  time  at  which  the  centerline  velocity  first  achieves  a  selected  percent 


of  Its  steady  state  value,  then  we  can  Identify  different  stages  In  the  propaga¬ 
tion  that  are  controlled  by  different  processes  and  therefore  are  likely  to 
scale  In  different  ways. 

For  distances  up  to  several  tens  of  nozzle  diameters  there  Is  likely  to  be 
a  compact  liquid  core  that  moves  at  the  Injection  velocity.  Its  most  forward 
surface,  however,  encounters  the  resistance  of  the  ambient  gas.  A  stagnation 
point  can  be  conceived  with  gas  forcing  the  liquid  to  flow  radially  away  from  It 
and  then  backward.  Thus  the  observed  tip  velocity  will  be  a  function  of  the 
Injection  velocity,  the  density  ratio,  and  the  actual  structure  of  the  core. 

Eventually  the  compact  core  disappears  and  ligaments  and/or  drops  encounter 
the  ambient  gas  first.  Since  drops  and  gases  exchange  forces  differently  than 
liquid  columns  and  gases,  the  apparent  tip  velocity  need  not  scale  as  It  did 
earlier.  At  sufficient  distance  from  the  nozzle,  most  of  the  Jet  momentum  Is 
with  the  entralni^d  gas  and  a  gaseous  head  vortex  moves  Into  the  gaseous  environ¬ 
ment  and  exchanges  forces  with  It  In  a  manner  which  Is  again  different  from  that 
of  drops  and  liquid  columns.  Thus  the  scaling  Is  likely  to  be  different  at  dif¬ 
ferent  stages  of  the  tip  propagation. 

Reitz  filmed  the  Initial  propagation  of  atomizing  Jets  but  he  was  con- 
cerned  with  an  overall  description  of  the  process  and  a  broad  set  of  conditions 
and  his  data  are  insufficient  to  determine  accurate  correlations  for  the  Initial 
tip  speed.  Far  from  the  Injector,  the  tip  velocity  is  about  70?;  of  the  velocitv 
given  by  Eqs.  11  and  5  (with  xg  neglected  with  respecc  to  x  or  estimated  usln. 
the  delta  of  Wu  et  al  There  Is  no  dependable  information  for  the  scaM 

of  the  intermediate  region  and  its  merging  with  the  other  two. 
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COMPUTATIONS 

.  A  good  test  of  the  degree  of  our  understanding  of  the  structure  of  these 
sprays  is  the  extent  to  which  we  can  predict  their  details* 

Navler-Stokes  equations  for  two-phase  flows  are  the  appropriate  conser¬ 
vation  equations  but  cannot  be  solved  with  adequate  resolution  for  the 

entire  flowfield  either  analytically  or  numerically  (3^).  Local  averaging  and 
the  adoption  of  semienplrlcal  equations  to  represent  the  effects  of  the 
neglected  details  are  necessary.  Even  so  no  attempt  has  been  made  to  compute 
the  flow  within  the  nozzle  to  predict  the  Initial  perturbations  required  by  the 
^^aerodynamic  theory  of  breakup  to  account  for  the  effects  of  the  nozzle  geometry. 

Many  difficulties  remain  even  if  one  considers  starting  the  computations  at 
the  nozzle  exit  from  some  arbitrary  initial  perturbations  of  the  liquid  gas 
interface*  Not  predicted  by  the  linear  stability  theory,  and  therefore  unknown, 
are  the  size  and  size  distribution  of  the  unstable  growths  at  breakup  and  the 
time  between  successive  ruptures*  Away  from  the  nozzle  exit  plane,  as  the 
generating  surface  moves  closer  to  the  axis  of  the  Jet,  there  are  questions  as 
to  what  gas  field  is  seen  by  the  liquid  Interface*  The  velocity  of  the 


entrained  gas  is  closer  to  that  of  the  generating  surface*  Thus  it  would  appear 


that  the  breakup  process  should  become  coupled  to  the  structure  of  the  two-phase 
flow  field  that  exists  between  the  presumed  intact  core  and  the  unperturbed 
outer  gas.  As  the  relative  velocity  between  liquid  and  gas  decreases  Inside  the 
Jet,  larger  drops  or  ligaments  or  blobs  should  be  formed;  Just  as  different 
breakup  modes,  and  larger  drops  are  found  when  the  injection  velocity  is 
decreased,  as  in  going  from  Fig.  2d  to  Fig.  2a.  Also,  coalescence  of  the  liquid 
fragments  can  be  expected  where  formation  occurs  due  to  the  locally  large  value 
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of  the  liquid  volume  fraction.  The  net  outcome  may  be  the  small  difference  be¬ 
tween  large  formation  and  coalescence  rates.  Thus  It  Is  clear  that  Eqs.  13-17 
may  provide  some  Information  about  the  outcome  of  the  breakup,  but  In  no  way  do 
they  give  all  that  Is  necessary. 

Faced  with  so  many  unknowns,  even  the  most  advanced  of  current  spray  com¬ 
putations  (3^),  are  Initiated  at  some  distance  from  the  nozzle  from  selected 
drop  size  and  velocity  distributions.  Then  Eqs.  15-17  are  used  to  give  limits, 
average  values,  and  functional  dependences.  The  results  should  become  Indepen¬ 
dent  of  the  details  of  the  Initial  conditions  at  large  distances  from  the 
nozzle.  This  tends  to  be  the  case,  for  example,  for  the  computed  size  of  the  ^ 
drops  when  drop  collisions  and  coalescence  are  Included  because  a  higher  rate  of 
coalescence  Is  found  from  smaller  Initial  drops  than  from  larger  ones.  Since, 
as  previously  discussed,  the  far  field  is  dominated  by  the  entrained  gas,  the 
only  necessary  conditions  for  the  accurate  computation  of  far  field  mean  quan¬ 
tities  are  knowledge  and  conservation  of  axial  momentum  and  proper  gas  phase 
turbulent  dlffuslvlty  (35),  Even  those  minimal  requirements  are  not  met 
without  some  care.  In  general  it  Is  difficult  accurately  to  compute  the  Injec¬ 
tion  momentum  of  an  actual  spray  even  when  Its  mass  flow  rate  Is  known,  because  ^ 
the  Injection  velocity  profile  Is  not  known  (Cm  In  E-;.  1  Is  not  know'n  even  If  Cq 
In  Eq.  2  Is  known),  k-e  models  can  be  tuned  to  give  the  far  field  dlffuslvlty 
of  incompressible  jets  If  such  dlffuslvlty  were  known  accurately.  But  this  is 
not  the  case  as  shown,  for  example,  by  varying  v-  chin  ±10%  when  different 
sets  of  experimental  data  are  used  (^3),  Succcssf-.  '.  computations  nf  fluc¬ 
tuations,  such  as  those  of  the  drop  velocity,  requii  that  the  model  be  accurate 
also  at  smaller  scales.  The  comparisons  (33)  g  9  adequate 
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reproduction  of  mean  axial  velocities  but  the  amplitude  of  the  drop  axial  velo¬ 
city  fluctuation  la  underestimated  by  a  factor  of  two.  Several  reasons  could  be 
advanced  for  the  specific  disagreement  of  Fig.  9  but,  ultimately,  current  quan¬ 
titative  knowledge  of  the  coupling  between  drop  motion  and  turbulence  gas  eddies 
Is  Incomplete  ^  broad  review  of  the  many  approaches  used  In  modeling 

sprays  and  other  two-phase  flows  Is  given  by  Faeth  (1^). 

SUMMARY 

What  Is  known  about  the  structure  of  non  vaporizing,  non  reactive  sprays 
from  single  hole  cylindrical  nozzles  can  be  summarized  with  the  help  of  Fig.  10. 

At  sufficiently  high  Injection  velocity  the  jet  Is  found  to  diverge  Imme¬ 
diately  at  the  nozzle  exit  where  fine  drops  are  observed.  The  outer  surface  of 
the  liquid  Is  disrupted  by  the  Interaction  with  the  ambient  gas  that  leads  to 
the  rapid  and  selective  growth  of  surface  waves  whose  Initial  amplitudes  are 
controlled  by  events  that  occur  within  the  nozzle.  This  view  allows  one  to  pre¬ 
dict  the  Initial  angle  of  the  spray.  The  parallel  prediction  of  the  size  of  the 
drops  thus  formed  also  compares  favorably  with  measurements  If  collisions  and 
coalesce  of  the  drops  after  their  formation  are  Included. 

The  disruption  of  the  liquid  column  eventually  reaches  the  axis  of  the  jet, 
i.e.  Is  complete,  because  only  Isolated  drops  are  found  downstream.  However  the 
geometry,  structure,  and  mode  of  breakup  of  the  core  of  the  Jet  are  not  known, 
but  there  are  Indications  that  the  length  of  the  Intact  core,  xi,  may  approach 
100  nozzle  diameters  (versus  10,  for  Incompressible  jets). 

While  the  Internal  breakup  continues,  gas  from  the  environment  Is  entrained 
rapidly,  the  entrained  gas  eventually  becomes  dominant  and  achieves  the  struc¬ 
ture  of  fully  developed  Incompressible  jets.  Within  It,  drops  tend  to  reach 
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equllibrluQ  buC  selectivity  since  both  drop  and  eddy  sizes  are  distributed.  The 
fully  developed  equilibrium  distance,  X2,  is  several  hundred  nozzle  diameters 
(versus  several  tens,  for  incompressible  Jets).  The  virtual  origin  of  this 
fully  developed  far  field,  xq.  Is  of  the  order  of  tens  of  nozzle  diameters 
(versus  order  of  nozzle  diameters,  for  incompressible  Jets). 

The  precise  relationship  between  xq,  xx,  and  X2  and  shape  of  the 
corresponding  boundaries  is  not  known  because  the  structure  of  the  transition 
region  is  not  known.  But  the  three  lengths  Increase  with  increasing  liquid-gas 
density  ratio. 

The  propagation  velocity  of  the  tip,  or  head,  of  these  sprays  is  a  fraction 
of  their  steady  state  centerline  velocity.  In  the  fully  developed  equilibrium 
region  this  fraction  is  about  70Z  but  in  the  development  region  it  is  not  known. 
Since  the  steady  state  centerline  velocity  scales  in  an  established  way  in  the 
fully  developed  equilibrium  region,  so  does  the  tip  velocity.  Experimental 
information  suggests  that  the  steady  state  centerline  velocity  and  the  tip  velo¬ 
city  scale  also  in  the  development  region,  at  least  approximately.  But  these 
Scaling  functions  are  determined  by  the  structure  of  the  development  region  and 
therefore  are  not  known  precisely. 
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Table  1:  Spray  Conditions  for  the  Drop  Velocity  Measurements 
of  Figs.  3,  8,  and  9  (from  Ref.  15) 


Series  P 


Pg  Pg/P£ 

Ap 

'^inj 

Nozzle 

X/d 

(MPa) 

(MPa) 

(m/s) 

d(wm)-L/d 

1.48  0.0256 

11.0 

127 

127-4 

600,800 

4.24  0.0732 

11.0 

127 

127-4 

300,400,500,600 

4.24  0.0732 

26.2 

194 

127-4 

400,500,600 

4.24  0.0732 

11.0 

149 

76-4 

300,600,700,800 

1.48  0.0256 

11.0 

125 

76-1 

300 

Liquid  ;  n-hexane, 

P£  -  665  kg/m^,  p£ 

-  3.2  X  10^ 

IJ.  s/m^. 

0£  -  1.84  X  102  N/ 

m.  Gas 

:  nitrogen. 

Room  temperature. 
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Table  2:  Spray  Conditions  for  the  Drop  Size  Measurements 
of  Figs.  6  and  7  (From  Refs.  26,  27) 


Series 

Nozzle 

L-d 

(pm) 

**8 

(MPa) 

Pg 

(kg/m3) 

Liquid 

Ap 

(MPa) 

<1 

0 

Cd 

Rej 

xlO”^ 

Wej 

xlO-"* 

1 

335-4 

1.48 

17.0 

n-C6H14 

1.38 

59.4 

0.92 

4.14 

4.26 

2 

335-4 

2.86 

33.0 

n-C6H14 

1.38 

52.4 

0.81 

3.65 

3.32 

3 

335-4 

1.48 

17.0 

n-C6H14 

3.45 

79.2 

0.78 

5.51 

7.58 

4 

335-4 

2.86 

33.0 

n-C6H14 

3.45 

92.8 

0.91 

6.46 

10.41 

5 

335-4 

1.48 

17.0 

n-C6H14 

6.90 

99.0 

0.69 

6.89 

11.85 

6 

335-4 

2.86 

33.0 

n-C6H14 

6.90 

111.0 

0.77 

7.73 

14.89 

7 

335-10 

1.48 

17.0 

n-C6H14 

3.45 

79.2 

0.78 

5.51 

7.58 

8 

127-4 

1.48 

17.0 

n-C6H14 

3.45 

78.1 

0.77 

2.06 

2.79 

9 

335-4 

1.48 

17.0 

n-C14H30 

3.45 

81.2 

0.86 

0.95 

6.31 

No  ol 
Dropa 

119 

109 

116 

111 

107 

117 

119 

114 

103 
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Development  of:  a)  a  turbulent  incompressible  Jet 
(water  in  water,  Reo  *  4  lO^);  and  b)  an  atomizing  spray 
(50*  glycerol  +  50%  water  into  N2) 

Some  regimes  of  breakup  of  liquid  jets 
Mean  axial  drop  velocity  and  amplitude,  skewness,  and  flat¬ 
ness  of  the  fluctuation  of  the  axial  drop  velocity  for  the 
atomizing  sprays  of  Table  1. 

Initial  angle  of  atomizing  Jets  versus  density  ratio  with 
fixed  nozzle  geometry.  Injections  of  glycerol-water,  water, 
hexane,  tetradecane  into  N2»  Hei  Xe»  Ar  at  liquid  pressures 
of  500-13,300  psi  and  D  -  254,  343,  and  610  pm.  Room  temp¬ 
erature  (20), 

Initial  angle  of  atomizing  Jets  versus  nozzle  length-to- 
diameter  ratio:  q  from  Ref.  20;  X  from  Ref.  19. 

Edge  of  spray  and  droplets  in  the  immediate  vicinity  of  the 
nozzle  exit  for  an  atomizing  Jet  (26-27), 

Average  diameter  of  drops  in  the  immediate  vicinity  of  the 
nozzle  exit  for  the  atomizing  Jets  of  Table  2:  O'neasured; 

#  expected  trends  from  the  supplemented  aerodynamic 
theory  (26-27), 

Measured^®)  mean  axial  drop  velocity  and  computed  (35) 
mean  axial  drop  and  gas  velocities  for  Spray  A  of  Table  1. 
Measured  (1®)  amplitude  of  the  fluctuation  of  the  axial  drop 
velocity  and  computed  (35)  amplitudes  of  the  fluctuations  of 
the  axial  drop  and  gas  velocities  for  Spray  C  of  Table  1. 
Schematic  structure  of  atomizing  Jets. 
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